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SECTION  I 

INTRODUCTION  AND  SUMMARY 


The  holographic  Head-Up  Display  (HUD)  offers  performance  improve¬ 
ments  over  conventional  HUD  systems  in  two  respects:  the  light 
losses  in  the  combiner  are  greatly  reduced,  and  the  f ield-of-view  is 
increased.  Asymmetries  that  result  from  the  off-axis  nature  of  the 
holographic  combiner,  however,  aggravate  image  abberations  and  cause 
a  significant  amount  of  geometric  distortion  in  the  imagery.  Al¬ 
though  it  is  possible  to  provide  distortion  correction  electronically 
at  the  display  source,  optical  correction  is  preferable  to  reduce 
the  complexity  of  the  required  electronics,  in  particular  where  the 
display  source  is  a  matrix  type  of  display  device.  In  previous  holo¬ 
graphic  HUD  systems,  much  of  the  distortion  introduced  by  the  holo¬ 
gram  has  been  corrected  by  the  relay  optics.  In  this  report,  we 
present  the  results  of  a  design  investigation  of  a  low-distortion 
holographic  HUD  optical  system  in  which  distortion  was  corrected  in 
the  holographic  combiner  as  well  as  in  the  relay  lens.  To  accomplish 
this,  we  assumed  the  combiner  to  be  a  holographic,  or  diffractive, 
optical  element  formed  with  analytically  defined  construction 
wavefronts . 

Analytically  defined  construction  wavefronts  offer  superior  flex 
ibility  and  ease  of  design  compared  to  previous  methods  of  design 
using  complicated  construction  beams.  In  the  new  method,  the  pnase 
of  each  construction  wavefront  at  the  hologram  surface  is  established 
directly  by  means  of  a  mathematical  description,  such  as  a  power 
series  expansion;  the  actual  construction  beam  is  not  designed  until 
the  hologram  design  is  complete.  The  analytical  description  of  the 
wavefronts,  furthermore,  permits  the  use  of  wavefronts  that  are  dif¬ 
ficult  or  impossible  to  realize  with  conventional  optical  systems. 
These  wavefronts  are  implemented  through  the  use  of  computer  gener¬ 
ated  holograms  in  the  construction  beams,  with  conventional  optical 
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elements  providing  the  major  wavefront  curvature  and  the  computer 
generated  hologram  adding  smaller  phase  perturbations  to  the 
wavefront. 

This  report  describes  a  research  program  that  included  an  inves¬ 
tigation  of  aspheric  holograms  made  with  analytically  defined  con¬ 
struction  beams  and  the  design  of  a  low-distortion  holographic  HUD. 
The  aspheric  hologram  investigation  developed  a  general  method  for 
designing  the  construction  wavefronts  of  a  holographic  optical  ele¬ 
ment  based  on  minimizing  the  mean-squared  phase  error  in  the  image 
formed  by  the  hologram  while  satisfying  the  Bragg  condition.  Consid¬ 
eration  was  given  to  the  physical  realizability  of  the  wavefronts 
found  by  the  design  method.  This  investigation  is  described  in  Sec¬ 
tion  [[  of  the  report. 

A  low-distortion  holographic  HUD  optical  system  was  designed  to 
fit  the  configuration  of  the  F-16  cockpit.  The  combiner  was  designed 
by  describing  the  phase  of  each  construction  wavefront  in  terms  of  a 
series  of  Legendre  polynomials.  The  optical  system  included  a  seven- 
element  relay  lens  to  form  an  intermediate  image  of  the  display, 
aberrated  in  such  a  manner  as  to  compensate  for  the  combiner  aberra¬ 
tions.  Reduction  of  both  distortion  and  image  errors  was  carried 
out  through  computer  optimization  of  the  hologram  and  relay  lens 
parameters.  The  design  procedures  are  described  in  Sections  III  and 
IV,  which  include  the  complete  design  parameters  of  the  optical 
system. 

Section  V  is  a  discussion  of  the  system  performance,  including 
distortion,  ray  errors,  and  hologram  efficiency.  Over  most  of  the 
instantaneous  f i eld-of-v iew,  the  distortion  met  the  goal  of  0.5  per¬ 
cent  or  less.  Ray  errors  tended  to  be  somewhat  higher  than  the 
goals,  with  accuracy  errors  generally  smaller  than  parallax  errors. 
The  hologram  efficiency  was  satisfactory  over  most  of  the  field-of- 
view  and  at  most  head  locations. 
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Because  the  construction  wavefronts  of  the  holographic  combiner 
were  described  as  mathematical  expressions  at  the  hologram  surface 
only,  in  Section  VI,  we  discuss  realization  of  such  construction 
beams,  including  a  discussion  of  the  requirements  on  the  computer 
generated  hologram. 

In  Section  VII,  we  present  our  conclusions  and  recommendations. 
A  major  conclusion  is  that  distortion  correction  is  feasible  in  a 
holographic  HUD  that  uses  an  aspheric  hologram  as  the  combiner.  The 
positive  results  of  this  study  suggest  continued  investigation  of 
the  aspheric  holographic  element  both  at  a  basic  level  and  as  applied 
in  optical  systems  such  as  the  head-up  display. 
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SECTION  II 

INVESTIGATION  OF  ASPHERIC  HOLOGRAMS 


This  section  summarizes  the  results  of  an  investigation  of  design 
methods  for  aspheric  holograms.  In  a  previous  research  effort  [1], 
it  was  shown  experimentally  that  a  holographic  lens  having  one  of 

its  recording  wavefronts  defined  by  a  general  power  series  expansion 
could  achieve  significantly  improved  performance  over  a  finite  field 
as  compared  to  a  conventional  holographic  optical  element  (HOE).  A 
conventional  HOE  is  defined  to  be  a  HOE  for  which  both  recording 
wavefronts  are  spherical,  as  depicted  by  Figure  la.  It  should  be 
noted  that  a  plane  wavefront  is  considered  to  be  a  spherical  wave- 
front  with  an  infinite  radius  of  curvature.  An  aspheric  HOE,  on  the 
other  hand,  is  defined  to  be  a  HOE  for  which  one  or  both  recording 
beams  are  nonspher  ical ,  as  shown  in  Figure  lb.  We  emphasize  that 

the  spherical  and  aspheric  descriptors  do  not  refer  to  the  substrate 
of  the  hologram,  which  in  either  case,  may  be  flat  or  curved. 

In  previous  holographic  HUD  investigations  [2],  aspheric  con¬ 

struction  wavefronts  have  been  obtained  by  the  introduction  of  con¬ 
ventional  optical  elements,  often  anamorphic,  into  the  construction 
beams.  In  this  report,  the  term  aspheric  HOEs  refers  to  elements 

made  with  construction  wavefronts  that  are  analytically  defined. 
For  these  wavefronts,  we  describe  the  phase  at  the  hologram  surface 
in  terms  of  a  mathematical  expression,  such  as  a  power  series 
expansion . 

Because  of  the  additional  degrees  of  freedom  available  in  the 
construction  beams,  aspheric  HOEs  represent  a  much  larger  class  of 
optical  elements  than  conventional  HOEs.  It  is  expected  that  with 
these  additional  degrees  of  freedom,  aspheric  HOEs  will  generally  be 
capable  of  improved  performance  in  imaging  applications  as  compared 
to  conventional  HOEs. 
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Object  Beam 


(a)  Conventional  Holographic  Optical  Element  With  a  Spherical  Object  Beam 


(b)  Aspheric  Holographic  Optical  Element  With  on  Aspheric  Object  Beam 

Figure  1.  Construction  Configuration  of  Conventional  and  Aspheric  Holographic 
Optical  Elements 


1.  GEOMETRIC  CONSIDERATIONS 


Consider  for  the  moment  a  HOE  which  is  recorded  as  shown  in 
Figure  2a  and  which  is  read  out  as  in  Figure  2b.  Here  we  have 
assumed  that  the  hologram  itself  is  the  entrance  pupil  of  the  system. 
In  addition,  we  assume  that  the  element  is  operated  over  a  collection 
of  field  angles  (i.e.,  the  reconstruction  wavefront  impinges  upon 
the  hologram  over  a  range  of  angles  and  correspondingly  produces 
image  wavefronts  over  a  range  of  angles).  J.R.  Fienup  has  shown  [3] 
that  the  hologram  phase  $(x)  which  minimizes  the  weighted  mean- 
squared  phase  error  is 


where  A  is  the  ideal  function  at  field  angle  n  and  W  is  its 

n  3  n 

corresponding  weighting  function.  The  hologram  phase  is  the  differ¬ 
ence  of  the  construction  wavefront  phases 

tf(x)  =  0Q(x)  -  0R(x)  (2) 

where  #Q(x)  is  the  object  wavefront  phase  and  tf>R(x)  is  the  ref¬ 
erence  wavefront  phase.  Conceivably,  the  reference  wavefront, 
0R(x),  could  be  chosen  to  be  a  wavefront  which  is  readily  gener¬ 

ated,  such  as  a  spherical  wavefront,  and  then  the  object  wavefront, 
A0(x),  could  be  computed  using  Eq.  (2)  to  produce  the  desired 

hologram  phase  A(x).  It  should  be  noted,  however,  that  the  choice 
of  ^r(x)  is  not  totally  arbitrary.  In  general,  any  realizable 

phase  function  must  obey  the  following  inequality, 

(3) 

in  order  that  it  may  represent  a  propagating  wavefront  of  wavelength 
k.  Violation  of  this  inequality  implies  that  the  wavefront  is 


(a)  Construction 


(b)  Reconstruction 


Figure  2.  Construction  and  Reconstruction  Configuration  of  a  Conventional  Reflection 
Holographic  Element 
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evanescent.  Both  dQ(x)  and  d^(x)  must  satisfy  inequality  (3). 
This  implies  that 


^0 

ad  3*R 

ax 

ax  ax 

(4) 


must  also  be  satisfied.  Oepending  upon  the  values  of  ad/ax,  this 
may  constrain  the  values  of  ad^/ax  that  the  reference  wavefront, 
dR(x),  may  have. 

The  direction  cosines  which  describe  the  direction  of  propagation 
of  a  wavefront  ray  may  be  derived  from  a  knowledge  of  the  phase  of 
the  wavefront  on  a  plane  surface  as  follows: 


x  ad 
1  "  2w  ax 


(5) 


m 


_  x_  ad. 
=  2ir  ay 


(6) 


n  =  ±(1 


2\M  2 
-  m  ) 


(7) 


where  £,  m,  n  are  the  x,  y,  z  direction  cosines,  respectively. 

In  terms  of  Eqs.  (1)  and  (2),  the  hologram  phase  has  an  x  direc¬ 
tion  cosine 


or 


£ 


x  ad 
2 n  ax 


I 


ad 

W  —  (x) 
n  ax  v  ' 


(8) 
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1  N 


adR 

ax 


(9) 


In  other  words,  if  we  think  of  the  hologram  phase  as  representing  a 
propagating  wavefront,  we  see  that  its  x  direction  cosine  is  just 
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the  weighted  average  of  the  difference  of  the  direction  cosines  of 
the  recording  wavefronts. 

2.  ENTRANCE  PUPIL  IN  FRONT  OF  HOE 

Next,  we  consider  readout  in  the  case  for  which  the  entrance 
pupil  is  in  front  of  the  HOE  as  depicted  in  Figure  3.  Intuition 
would  tell  us  that  this  geometry  is  capable  of  greater  performance 
over  extended  fields-of-view  because  at  the  extreme  field  angles, 
the  readout  wavefront  illuminates  totally  different  sections  of  the 
hologram.  By  appropr iately  recording  the  hologram  in  these  different 
sections,  the  readout  geometry  may  more  closely  correspond  to  the 
recording  geometry  and  thus  better  performance  may  be  achieved.  An¬ 
other  way  of  looking  at  this  is  as  follows.  At  each  position  x  on 
the  hologram,  a  phase  equal  to  the  weighted  average  of  the  difference 
of  the  recording  wavefront  phases  is  recorded.  The  RMS  error  of  the 
average  phase  will  be  proportional  to  the  ray  direction  errors  (aber¬ 
rations)  upon  readout  at  various  field  angles.  In  general,  the 
larger  the  range  of  field  angles  which  illuminate  hologram  position 
x,  the  larger  will  be  the  RMS  error.  For  the  case  of  the  stop  lo¬ 
cated  at  any  position  other  than  at  the  hologram  plane,  the  range  of 
field  angles  which  illuminate  hologram  position  x  will  be  less  than 
the  total  field  of  view  of  the  system  and  thus  the  RMS  error  will  be 
reduced. 

Let  us  now  state  what  this  means  in  terms  of  Eq.  (1).  We  still 
want  Eq.  (1)  to  be  satisfied  at  each  hologram  position  x;  now,  how¬ 
ever,  the  range  of  field  angles  over  which  the  summation  is  computed 
is  a  function  of  the  hologram  position,  x.  This  field  angle  depend¬ 
ence  on  x  can  be  conveniently  incorporated  into  the  weighting 
function . 

Assume  the  readout  geometry  of  Figure  3  with  an  entrance  pupil 
of  extent  Xp  oriented  parallel  to  the  hologram  and  illuminated  by 


plane  waves  over  an  extended  field  of  view.  Consider  a  ray  beginning 
at  the  center  of  the  entrance  pupil  and  intercepting  the  hologram  at 
Xp  This  defines  a  field  angle  in  terms  of  its  hologram  inter¬ 
cept.  For  each  field  angle,  we  can  define  a  hologram  phase  as  a 
function  of  hologram  coordinate.  This  corresponds  to  the  difference 
between  the  phases  of  the  readout  wavefront  at  the  specified  field 
angle  and  its  corresponding  desired  image  wavefront.  Thus,  the  holo¬ 
gram  phase  may  be  described  as  a  two-dimensional  function  of  field 
angle  (as  specified  by  hologram  intercept)  and  hologram  coordinate, 

<bh  =  x)  (10) 

where  and  x  are  the  field  angle  intercept  and  hologram  coordi¬ 
nate,  respectively.  If  we  now  also  define  a  weighting  function, 
W(x),  such  that 

W(x)  =  0  for  j  x  {  >  ( 11 ) 

and 

oo 

f  W(x)  dx  =  1  (12) 


then  we  can  define  the  composite  hologram  phase  as 

00 

d(x)  =  /  W(x  -  xf)dh(xf,  x)  dxf  (13) 

CO 

Once  the  composite  hologram  phase  of  Eq.  (13)  is  determined,  per¬ 
haps  numerically,  the  task  remains  as  before  to  choose  object  and 
reference  wavefronts  to  realize  this  phase. 

The  analysis  we  have  performed  to  this  point  has  assumed  a  thin 
hologram  structure.  In  particular,  the  efficiency  of  the  diffracted 
wavefront  is  not  considered.  In  light  of  the  fact  that  most  practi¬ 
cal  applications  will  involve  thick  holograms,  Bragg  diffraction 
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effects  must  be  considered.  The  following  section  considers  an  al¬ 
ternate  approach  to  the  diffraction  analysis  which  includes  Bragg 
diffraction  requirements. 


3.  VECTOR  ANALYSIS  OF  HOLOGRAPHIC  IMAGING 

An  alternate  approach  to  analyzing  the  diffraction  of  a  holo¬ 
graphic  grating  is  through  a  vector  analysis  as  described  by  Spencer 
and  Murty  [4]  and  by  Welford  [5].  This  technique  is  used  extensively 
in  hologram  ray  tracing  programs,  including  HOAD. 

Let  the  reconstruction  and  image  wavefronts  be  described  at  any 

position  by  rays  perpendicular  to  the  wavefronts  in  the  direction  of 

propagation.  Let  the  rays,  in  turn,  be  described  by  unit  vectors,  C 
■+ 

and  I  in  a  cartesian  coordinate  space.  The  components  of  the  unit 
vectors  are  simply  the  direction  cosines,  i.e.. 


and 


» 

C  =  ix 
c 


+  V  +  V 


(14) 


^  A  A 

I  =  +  mjy  +  n jZ  (15) 

A  A  A 

where  x,  y  and  z  are  unit  vectors  forming  a  right-handed  coordinate 

^  ^  a 

system  such  that  x  x  y  =  z.  In  addition,  let  the  grating  be 
described  by  its  surface  period  (fringe  spacing),  d  and  orienta- 
tion,  q,  where  q  is  a  unit  vector  parallel  to  the  fringes  that  lies 
:u  the  plane  of  the  surface.  Then  the  following  relationship  holds 

(I  -  C)  x  n  =  (xc/ds)q  (16) 

A 

where  n  is  a  unit  vector  normal  to  the  surface  and  x  is  the  re- 

c 

construction  free-space  wavelength. 

Let  us  now  assume  that  the  grating  was  recorded  interferometri- 
cally  by  two  wavefronts  described  by  vectors  0  and  R.  The  recording 
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grating  will  have  fringe  planes  with  orientation  and  separation  as 
follows, 


>  .  »  » 

G  =  0  -  R  (17) 


d 


G 


(18) 


where  G  is  a  vector  normal  to  the  fringe  plane,  dg  is  the  fringe 
plane  spacing,  and  \Q  is  the  construction  free-space  wavelength. 
The  grating  vector  component  which  is  parallel  to  the  fringes  and 
lies  in  the  surface  tangent  plane  is  given  by 


Q  =  (0  -  R )  x  n 

and  thus,  the  corresponding  surface  grating  spacing  is 


d  = 


5  151  1(0  -R)  x  nl 

Combining  Eqs.  (16)  and  (19),  we  find 

(I  -  C)  x  n  =  ^  |  (0  -  R)  x  n|q 


but 


»  * 

a  (0  -  R)  x  n 

ho -b)  « ;i 

so 

»  >  X  >  > 

(I  -  C)  x  n  =  (0  -  R)  x  n 

xo 


(19) 


(20) 


(21) 


(22) 


This  is  the  standard  vector  ray  tracing  equation  for  holographic 
gratings  used  by  many  ray  tracing  design  programs. 
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4.  BRAGG  EFFICIENCY  CONSIDERATIONS 


Equation  (22)  describes  the  geometric  relationship  between  the 
holographic  recording  and  readout  wavefronts.  The  amplitude  of  the 
image  wavefront,  however,  is  a  complex  function  of  the  recording 
geometry,  readout  geometry  and  material  properties  [6].  It  is  suf¬ 
ficient  for  the  purposes  of  this  analysis,  however,  to  note  that  for 
peak  efficiency,  the  Bragg  condition  is  satisfied  when  the  recording 
grating  structure  is  identical  in  fringe  spacing  and  orientation  to 
that  which  would  be  recorded  by  the  desired  readout  wavefronts.  In 
equation  form,  we  require  that 


and 


x„  x„ 

o  c 

6  =  IT  -  Cl 


|0  -  itl 


0-R  _  ~  _  I  -  C 


lo-Rl 


IT  -£I 


(23) 


(24) 


where  dg  is  the  fringe  spacing  and  g  is  a  unit  vector  perpendicu¬ 
lar  to  the  fringe  planes.  The  above  two  equations  define  orthogonal 
measures  and  so  we  may  combine  them  into  a  single  equation 


g  0  -  R  lO  -  Rl  I  -  C  ll  -  Cl 

vro-v  xo  ‘if -ci'  xc 

or  finally 


(25) 


>  *  x  >  » 

0  -  R  =  -2-  (I  -  C)  (26) 

xc 

This  is  the  vector  Bragg  constraint  which  must  be  satisfied  to 
achieve  maximum  diffraction  efficiency  with  the  desired  imaging 
properties. 
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It  is  interesting  to  compare  the  imaging  Eq.  (22)  with  the  com¬ 
bination  Bragg  and  imaging  Eq.  (26)  just  developed.  Equation  (22) 
is  identical  to  Eq.  (26)  except  that  both  sides  of  Eq.  (22)  are 
crossed  with  the  surface  normal,  n.  In  geometric  terms,  Eq.  (22) 
relates  the  wavefront  components  in  the  tangent  plane  of  the  surface, 
whereas  Eq.  (26)  relates  the  wavefronts  in  a  volumetric  sense.  Al¬ 
ternatively,  we  may  say  that  Eq.  (26)  considers  the  grating  structure 
within  the  volume  of  the  hologram,  whereas  Eq.  (22)  considers  the 
grating  structure  only  upon  the  surface  of  the  hologram.  It  follows 
algebraically  that  if  Eq.  (26)  is  satisfied,  then  Eq.  (22)  is  satis¬ 
fied  also. 

5.  REALIZABLE  WAVEFRONTS 

As  stated  previously,  any  two-dimensional  phase  function  on  a 
plane  tf(x,y),  describes  a  realizable  wavefront,  provided  the  phase 
does  not  vary  too  rapidly  in  comparison  to  the  propagation  constant 
(see  inequality  (3)).  It  also  follows  that  given  a  phase  function 
which  satisfies  inequality  (3),  the  ray  direction  cosines  at  any 
point  on  the  wavefront  may  be  computed  from  the  partial  derivatives 
of  the  phase  function  as  described  in  Eqs.  (5)  -  (7).  From  a  purely 
mathematical  standpoint,  we  require  that  the  function  0(x,y)  be  both 
continuous  and  differentiable.  Obviously,  two-dimensional  polynomial 
functions  satisfy  both  of  these  mathematical  constraints  and  are 
easily  differentiated  algebraically.  Other  functions,  such  as 
weighted  sums  of  phase  functions,  may  also  satisfy  the  continuity 
and  differentiability  constraints,  but  may  require  numerical  methods 
to  approximate  the  partial  derivatives.  In  this  manner,  analytically 
defined  wavefronts  can  be  described  as  phase  functions  and  ray  trac¬ 
ing  methods  may  be  used  (as  in  HOAD)  to  analyze  the  performance  of 
such  a  HOE  in  an  optical  system. 

We  now  consider  the  inverse  problem,  i.e.,  given  two  functions 
(each  two-dimensional)  which  represent  the  partial  derivatives  of  a 


15 


phase  function  defined  on  a  plane,  how  do  we  recover  the  phase  func¬ 
tion  and  perhaps,  more  importantly,  how  do  we  make  sure  that  the 
phase  function  is  realizable.  From  vector  integral  calculus  [7],  we 
find  that 


*(x,y) 


(27) 


provided  the  line  integral  is  independent  of  path.  In  other  words, 
if  the  two-dimensional  funtions,  ad/ax  and  ad/ay,  are  such  that  the 
integral  of  Eq.  (27)  has  the  same  value  regardless  of  the  path  of 
integration  from  (0,0)  to  (x,y),  then  d(x,y)  is  single  valued  and  in 
fact  atf/ax  and  afll/ay  are  the  partial  derivatives  of  tf(x,y).  Inde¬ 
pendence  of  path,  however,  is  not  a  property  which  can  be  numerically 
verified,  so  we  need  some  other  method.  Fortunately,  a  second 
theorem  of  vector  integral  calculus  provides  an  answer.  It  states 
that  if  atf/ax  and  atf/ay  have  continuous  derivatives  in  a  simply  con¬ 
nected  domain  and  if 


2  ? 

a  d  _  31  _  ann 

ayax  ~  axay  r  ay  =  ax 


(28) 


then  the  integral  of  Eq.  (27)  will  be  independent  of  path.  The  re¬ 
quirement  that  the  derivatives  be  continuous  is  a  reasonable  one 
since  the  derivatives  are  proportional  to  the  ray  direction  cosines 
which  must  vary  in  a  smooth  manner  for  any  reasonable  wavefront.  In 
general,  we  consider  only  simply  connected  domains  (i.e.,  we  do  not 
allow  holes  within  the  area  of  the  hologram). 

Thus,  if  the  direction  cosines  (or  atf/ax  and  30/ay)  are  chosen 
such  that  Eq.  (28)  is  satisfied,  then  Eq.  (27)  may  be  used  to  deter¬ 
mine  the  corresponding  two-dimensional  phase  function.  Ultimately, 
it  is  the  phase  function  that  we  want  in  order  to  synthesize  a  com¬ 
puter  generated  hologram  that  will  produce  the  desired  wavefront. 

At  this  point,  we  make  some  observations  which  may  be  useful 
later.  First,  if  we  have  two  realizable  phase  functions,  d^(x,y) 
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and  *2(x,y),  then  any  linear  combination  of  these  phase  functions 
is  also  a  realizable  phase  function,  because  the  integration  of  Eq. 

(27)  is  a  linear  operation.  For  example,  dj(x,y)  =  0-|(x,y)  - 
2$£(x,y)  is  a  realizable  phase  function.  In  a  similar  manner,  if 
we  have  two  sets  of  direction  cosines  which  satisfy  Eq.  (28),  then 
any  linear  combination  of  the  sets  will  also  satisfy  Eq.  (28),  be¬ 
cause  the  derivatives  in  Eq.  (28)  are  linear  operations.  This  means 
that  the  linear  combination  may  be  used  in  Eq.  (27)  to  arrive  at  a 
valid  combination  phase  function.  For  example, 

i,  =  J'1  -  2?2  m^  =  m^  -  2m,,  (29) 

will  satisfy  Eq.  (28)  provided  and  m2  satisfy  Eq. 

(28)  separately. 

The  final  observation  we  wish  to  make  is  that  general  two- 
dimensional  polynomials,  such  as. 


m  n 

x,y)  =  1C  aijx1yJ  (3°) 

i  =0  j=0 


are  realizable  phase  functions,  as  are  all  linear  combinations  of 
such  polynomials.  Furthermore,  any  polynomial  pairs  which  satisfy 
Eq.  (28)  also  can  be  used  with  Eq.  (27)  to  obtain  a  realizable  phase 
function.  For  example,  suppose  we  have 

t  =  2x  +  y  m  =  x  +  2y  (31) 


Then 


and 


a ;  |  am 

ay  “ax 


(32) 


fx,y 

tf(x, y)  =  ^7  [  tdx  +  m  dy] 

cJo,o  A 


(33) 
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6.  IMPLEMENTATION 

We  will  now  specify  a  method  whereby  Eq.  (26)  may  be  used  to 
determine  the  recording  wavefronts  to  obtain  the  desired  performance. 
We  first  break  Eq.  (26)  into  its  components. 


xQ 

~  *r  =  ( _  ^c  ^ 

c 

%  -  m»  ■  r  <mi  -  "c1 

c 


(34) 

(35) 


n0  -  nR  =  r  (ni  -  "c>  (36) 

c 

Presumably,  we  know  what  ?p  m^,  n^,  S.  c,  nc,  xQ  and  xc 
are  desired  to  be  (as  functions  of  the  hologram  coordinates,  x  and 
y).  If  we  are  judicious  and  choose  the  pairs  Up  and  ^c> 
mc  such  that  Eq.  (28)  is  satisfied  for  each  pair  (perhaps  by 
choosing  them  as  appropriate  polynomials),  then  we  know  that  the 
right-hand  sides  of  Eqs.  (34)  -  (36)  will  represent  a  realizable 
phase  function.  It  then  remains  to  choose  the  left-hand-side  func¬ 
tions  such  that  they  are  realizable  and  equal  to  the  right-hand  re¬ 
sult.  One  possible  method  would  be  to  choose  Jt  ^  and  mR  such 
that  they  represent  some  known  realizable  phase  function  (such  as  a 
polynomial)  and  then  by  the  linearity  we  spoke  of  earlier,  and 
mQ  will  automatically  represent  a  realizable  phase  function.  This 
method,  however,  removes  some  degrees  of  freedom  from  the  problem. 
The  result  may  be  that  the  recording  beams  ft  and  ft  may  not  be  the 
pair  with  the  lowest  space-bandwidth  product  that  could  be  used.  A 
second  possibility  is  to  let 


1  xo  ,  V 

&  0  ^R  =  2\  ^1  ^c 

c 


(37) 
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(38) 


"o  -  ”r  *  r  (7  (mi  -"c> 


1  xo  , 

no  -  nR  -  2  xc  (nI  "  V 


(39) 


■+  -> 

in  which  case  0  and  R  will  have  equal  space  bandwidth  products. 

Finally,  once  0  and  R  have  been  established,  it  is  necessary  to 

evaluate  Eq.  (27)  either  algebraically  if  0  and  R  are  described  by 

■+ 

an  algebraic  equation,  or  numerically  if  0  or  R  are  described  by 
samples . 


\ 
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SECTION  III 

INITIAL  DESIGN  INVESTIGATIONS 


The  objectives  of  the  initial  design  investigations  were  to  lo¬ 
cate  the  major  system  components  and  to  establish  an  aspheric  com¬ 
biner  design  suitable  as  a  starting  point  in  the  subsequent  system 
design . 

1.  INITIAL  DESIGN  OF  THE  HUD  SYSTEM 

The  objective  of  the  initial  design  was  to  locate  the  major  com¬ 
ponents  of  the  optical  system,  and  to  establish  the  first-order  con¬ 
figuration  of  the  construction  beams.  We  assumed  a  system  configura¬ 
tion  of  the  form  shown  in  Figure  4.  The  initial  design  was  based  on 
satisfying  three  conditions.  The  first  was  the  requirement  that  the 
optical  system  fit  into  the  F— 1 6  cockpit,  which  constrained  the  loca¬ 
tions  of  the  combiner  and  fold  mirror.  The  second  condition  was  that 
the  holographic  combiner  have  high  diffraction  efficiency,  which 
constrained  the  general  location  of  the  construction  beams.  The 
third  condition  on  the  initial  design  was  that  the  hologram  design 
control  astigmatism  to  a  first  order,  and  where  possible,  minimize 
other  image  aberrations. 

a.  SYSTEM  LAYOUT 

The  major  objective  of  the  system  layout  was  the  placement  of 
the  combiner  and  fold  mirror  in  the  cockpit.  This  task  was  governed 
by  the  requirements  of  the  optical  system  and  the  configuration  of 
the  cockpit.  The  optical  system  requirements  that  affected  the  sys¬ 
tem  layout  were  the  size  of  the  exit  pupil  (3.0  inches  vertical  by 
4.5  inches  horizontal),  and  the  instantaneous  and  vertical  fields-of- 
view,  which  we  defined  as  follows. 

The  instantaneous  f ield-of-view  (IFOV)  measures  12°V  by  20°H;  it 
is  measured  from  the  center  of  the  exit  pupil  and,  in  the  vertical 
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plane,  is  defined  by  the  horizontal  and  a  ray  at  an  angle  of  12°  be¬ 
low  it.  In  the  horizontal  direction,  the  IFOV  is  defined  by  rays  at 
angles  of  *10°  with  the  vertical  plane. 

The  total  f ield-of-view  (TFOV)  is  20°,  *  10°  about  the  center  of 
the  field-of-view.  In  the  F-16  cockpit,  the  TFOV  is  defined  in  the 
vertical  plane  by  a  ray  from  the  top  of  the  exit  pupil  at  an  angle 
of  16°  below  the  horizontal  and  a  ray  up  from  the  bottom  of  the  exit 
pupil  at  a  4°  angle.  The  TFOV  is  circularly  symmetric. 

We  also  define  the  optical  axis  to  pass  through  the  center  of 
the  exit  pupil  at  the  center  of  the  f ields-of-view;  it  therefore 
points  downward  6°  with  respect  to  the  horizontal,  toward  the  nose 
of  the  aircraft. 

The  cockpit  boundaries  which  affect  the  HUD  system,  shown  in 
Figure  5,  are  the  following:  the  windshield  clearance  line  which  is 
located  parallel  to  and  3/4"  inside  the  windshield;  the  over-the- 
nose-vision-line  (OWL)  which  extends  down  from  the  center  of  the 
exit  pupil  at  15°  from  the  horizontal;  and  the  ejection  line.  The 
combiner  could  not  be  placed  forward  of  the  windshield  clearance  line 
and  no  part  of  the  system  could  block  the  ONVL.  We  allowed  some  vio¬ 
lation  of  the  ejection  line  on  the  grounds  that  the  optical  module 
was  sufficiently  limited  in  lateral  dimension  so  as  to  allow  adequate 
clearance  on  each  side.  The  optical  module  as  shown  in  Figure  5 
determines  boundaries  for  the  fold  mirror,  relay  lens,  and  CRT. 

The  combiner  was  positioned  as  far  forward  in  the  cockpit  as  was 
consistent  with  both  the  vertical  field-of-view  requirements  and  the 
available  space  between  the  windscreen  clearance  line  and  the  ONVL. 
The  folding  mirror,  on  the  other  hand,  was  positioned  as  far  aft  as 
the  optical  module  would  allow  in  order  to  minimize  the  angle  between 
the  chief  ray  incident  at  the  center  of  the  hologram  before  and  after 
diffraction  by  the  hologram,  in  other  words,  to  minimize  the  off-axis 
nature  of  the  hologram. 
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The  hologram  was  tilted  such  that  the  chief  ray  incident  at  the 
center  of  the  hologram  was  diffracted  at  an  equal  and  opposite  angle 
of  incidence,  as  though  it  had  been  reflected  by  classical  reflec¬ 
tion.  This  geometry  minimized  chromatic  effects  resulting  from  the 
spectral  bandwidth  of  the  phosphor  and  tended  to  minimize  distortions 
arising  from  asymmetries  in  the  grating  geometry. 

The  size  of  the  hologram  was  determined  by  the  eye  relief,  the 
f ields-of-view,  the  size  of  the  exit  pupil,  and  to  some  extent,  the 
curvature  of  the  hologram.  The  size  of  the  folding  mirror  required 
to  accommodate  all  rays  also  depended  on  these  parameters  and  in 
addition  on  the  focal  length  of  the  hologram,  the  separation  distance 
between  hologram  and  mirror,  and  the  space  allowed  by  the  optical 
modu 1 e . 

The  tilts  of  the  hologram  and  fold  mirror  were  adjusted  so  that 
a  ray  from  the  center  of  the  exit  pupil  in  the  center  of  the  field- 
of-view  would,  after  reflection  by  the  fold  mirror,  be  in  the  approx¬ 
imate  center  of  the  optical  module  parallel  to  the  ONVL. 

The  relay  lens  was  located  110  mm  from  the  fold  mirror,  centered 
between  the  top  and  bottom  of  the  optical  module.  Location  of  the 
CRT  was  put  off  until  the  relay  lens  parameters  were  defined  and  un¬ 
til  the  intermediate  image  requirements  were  established. 

b.  FIRST  ORDER  COMBINER  DESIGN 

The  objective  of  the  first  order  combiner  design  was  a  conven¬ 
tional  holographic  combiner  design  that  we  could  use  as  a  starting 
place  for  the  aspheric  holographic  combiner  design.  Requirements  on 
the  first  order  design  were  a  combiner  having  high  diffraction  effi¬ 
ciency  over  the  entire  vertical  field-of-view  and  having  an  image 
with  first  order  astigmatism  corrected.  A  major  factor  in  the  com¬ 
biner  design  is  the  wavelength  shift:  we  assumed  a  construction 
wavelength  of  514.5  nm  and  a  HUD  CRT  (reconstruction)  wavelength  of 
544  nm. 
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The  initial  combiner  design  proceeded  in  two  steps.  In  the 
first,  we  assumed  spherical  construction  beams,  with  parameters 
chosen  to  give  a  reasonable  focal  length  and  to  satisfy  the  Bragg 
condition.  In  the  second  step,  we  added  cylindrical  lenses  to  the 
construction  beams  to  correct  first  order  astigmatism  introduced  by 
the  hologram. 

The  configuration  of  Figure  4  called  for  a  reflection  type  dif¬ 
fractive  element,  which  is  constructed  with  the  point  sources  of  both 
construction  beams  lying  on  the  same  side  of  the  combiner,  with  one 
beam  convergent  and  the  other  divergent.  The  requirement  that  the* 
central  ray  be  diffracted  at  an  angle  equal  and  opposite  to  its  in¬ 
cident  angle  placed  the  construction  beam  point  sources  on  radii  mak¬ 
ing  equal  and  opposite  angles  with  the  normal  to  the  hologram  at  its 
center  (the  hologram  center  was  defined  to  be  at  the  intersection  of 
the  optical  axis  with  the  combiner). 

The  radii  of  the  beams,  Rq  and  Rf,  are  related  to  the  focal 
length  of  a  flat  combiner  by  the  equation 


where  f  is  the  focal  length,  \Q  is  the  wavelength  used  during  re¬ 
construction  and  \Q  is  the  wavelength  used  during  recording.  The 
curvature  of  the  combiner  also  affected  its  focal  length  but  this 
was  later  compensated  for  by  scaling  both  radii  by  a  suitable  amount. 
The  choice  of  RQ  and  Rr  was  based  on  the  first  order  imaging 
properties  of  the  HUD  optical  system.  Since  the  relay  lens  and  the 
exit  pupil  are  conjugate  points  with  respect  to  the  hologram,  reason¬ 
able  locations  for  point  sources  used  to  form  construction  beams  for 
holograms  were  the  center  of  the  relay  lens  and  the  center  of  the 
exit  pupil.  Although  other  pairs  of  conjugate  points  could  have  been 
selected  to  locate  the  construction  beam  point  sources,  these  points 
were  preferable  because  they  are  optimum  for  maintaining  a  high  dif¬ 
fraction  efficiency  over  the  field-of-view. 
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The  wavelength  shift  from  514.5  nm  of  the  laser  to  544  nm  of  the 
P-43  phosphor  required  a  significant  modification  of  the  construction 
beams  in  order  to  satisfy  the  Bragg  condition.  At  the  center  of  the 
hologram,  we  calculated  the  construction  beam  angles  required  to  es¬ 
tablish  fringes  with  correct  spacing  such  that  on  reconstruction, 
the  ray  incident  at  the  center  of  the  hologram  from  the  center  of 
the  exit  pupil  would  satisfy  the  Bragg  condition.  The  construction 
beam  angles  were  given  by  solving 

sin  15—  =  ^ —  sin  ^  (41) 


where  9  is  the  angle  between  the  beams  within  the  hologram  medium, 
and  x  is  the  wavelength  (in  air).  If  the  unprimed  quantities  refer 
to  the  hologram  during  reconstruction  (i.e.,  in  the  HUD  optical  sys¬ 
tem),  then  the  primed  quantities  are  the  corresponding  parameters 
for  the  hologram  construction. 


In  our  design,  the  ratio 


x'  0.544 
x  =  0.5)45  " 


1.057 


For  o  =  30° , 
then  e'  =  42.765° 


Since  the  reference  and  object  wavefronts  were  spherical,  this 
analysis  was  correct  only  at  the  hologram  center;  at  other  locations 
on  the  hologram,  the  construction  rays  do  not  make  equal  and  opposite 
angles  with  respect  to  the  hologram  surface  normal.  The  analysis 
was  sufficient,  nevertheless,  to  locate  the  construction  beams  to  a 
first  order,  leaving  refinements  to  a  later  step. 

The  construction  beam  lengths  were  also  modified  to  accommodate 
the  wavelength  shift.  The  reference  beam  length  became  the  distance 
from  the  exit  pupil  to  the  hologram  multiplied  by  the  ratio  of 
reconstruction-to-construction  wavelengths.  Similarly,  the  object 
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beam  length  was  made  the  distance  from  the  hologram  to  the  relay  lens 
multiplied  by  the  wavelength  ratio. 


The  first  step  in  the  combiner  design  was  completed  by  setting 
Rq  =  300  mm  and  Rr  =  500  mm.  As  hologram  curvature  is  a  useful 
independent  design  variable,  we  chose  a  substrate  radius  of  550  mm 
for  the  initial  design.  These  parameters  resulted  in  a  hologram 
focal  length  in  the  vertical  plane  of  about  223  mm,  with  approxi¬ 
mately  55  mm  separating  the  vertical  and  horizontal  foci. 

The  second  step  of  the  initial  combiner  design  was  to  add  a 
cylinder  lens  to  each  construction  beam  to  correct  astigmatism  intro¬ 
duced  by  the  hologram.  As  the  initial  design  was  rather  preliminary, 
no  attempt  was  made  to  find  the  optimum  location  and  power  of  the 
cylinder  lenses  once  the  astigmatism  had  been  corrected  to  a  first 
order . 

Figure  6  shows  the  configuration  of  the  construction  beams  that 
resulted  from  the  initial  design  investigation.  One  beam  is  diverg¬ 
ing  through  the  construction  beam  optical  system  as  shown.  The  other 
beam  is  considered  to  be  converging,  passing  from  the  combiner 
through  its  cylinder  lens  to  a  point  focus.  We  design  the  latter 
beam  by  tracing  rays  from  the  point  source  to  the  hologram,  where 
the  ray  directions  are  reversed  by  the  raytrace  program.  This  design 
resulted  in  a  combiner  with  a  focal  length  of  193  mm,  and  which  is 
predicted  to  maintain  over  90  percent  diffraction  efficiency  over 
the  vertical  field. 

2.  ASPHERIC  COMBINER 

The  aspheric  combiner  design  was  based  on  the  analytical  defini¬ 
tion  of  the  phase  of  each  construction  beam  at  the  hologram.  The 
objective  of  the  combiner  design  was  an  aspheric  hologram  with  the 
correct  first  order  imaging  properties,  corrected  astigmatism,  and, 
if  possible,  low  coma  and  low  distortion.  The  procedure  began  by 
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replacing  the  conventionally-defined,  anamorphic  construction  beams 
discussed  in  the  previous  section  with  analytically  defined  construc¬ 
tion  beams  having  the  approximate  forms  of  the  anamorphic  wavefronts. 
We  then  manually  adjusted  the  parameters  defining  the  construction 
wavefronts  to  improve  the  hologram  performance,  and  finally  optimized 
the  construction  wavefronts  using  the  optimization  routines  of  the 
HOAD  program  [8].  As  it  turned  out,  the  manual  adjustment  of  the 
construction  wavelength  parameters  was  sufficient  to  establish  a 
hologram  design  with  which  to  begin  system  integration  and 
optimization . 


a.  ANALYTICALLY  DEFINEO  WAVEFRONTS 


In  all  of  our  investigations,  the  phase  0(x,y)  of  each  construc¬ 
tion  beam  at  the  hologram  was  defined  in  terms  of  an  expansion  in 
Legendre  polynomials,  expressed  as 

9  9 

«x,y)  ,  £  £  C(JL((*r)Lj(yr)  (42) 

i=0  j=0 


in  which  x  and  y  are  the  relative  coordinates 
r  Jr 


x  -2 x  _2y 
r  "  0  ’  yr  ~  TT 


where  x  and  y  correspond  to  vertical  and  horizontal  coordinates  of 
the  hologram,  respectively,  C,.  are  coefficient  weighting  values 

■  J 

specified  by  the  user,  and  D  is  the  hologram  diameter.  Legendre 
polynomials  are  attractive  because  they  are  orthogonal,  which  means 
that  changing  the  coefficient  of  one  polynomial  does  not  affect  the 
weighting  of  the  other  polynomials.  Some  examples  of  Legendre  poly¬ 
nomials  are  given  below. 

Lj(x)  =  x  (43a) 
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l2(x)  =  (3x2  -  1)  (43b) 

L3(x)  =  \  (5x3  -  3x)  (43c) 

L4(x)  =  ^  (35x4  -  30x2  +  3)  (43d) 

Notice  that  the  magnitude  of  the  polynomial  is  unity  when  x  =  ±1 , 


with  the  sign  of  the  polynomial  at  x  =  -1  depending  on  the  order  of 
the  polynomial.  It  is  because  of  this  property  that  we  used  relative 
hologram  coordinates  in  Eq.  (42). 

As  an  example  of  the  use  of  the  Legendre  polynomial  expansion  to 
describe  a  wavefront,  let  us  consider  the  case  of  a  spherical  wave- 

front.  We  assume  the  geometry  shown  in  Figure  7,  with  the  hologram 

in  the  x-y  plane,  illuminated  by  a  spherical  wave  diverging  from  the 
point  xQ,  0,  zQ  located  a  distance  R  from  the  center  of  the 

hologram.  The  phase  tf(x,y)  of  the  wavefront  at  the  hologram  to 

fourth  order  is  expressed  as 

tf(x,y)  =  C10L1  'xr^ 

+  C20L2(xr)  +  co2L2^yr^ 

+  C12Ll^xr^L2^yr^  +  C30L3^xr^ 

+  C40L4^Xr^  +  C22L2^xr^L2^yr^  +  C04L4^yr^  ^44^ 

where,  again, 


r  IT 

The  coefficients  are  given  by 


x  _  =  ^  ,  yr  =  / 
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b.  DESIGN  OF  THE  ASPHERIC  COMBINER 

The  aspheric  combiner  design  began  by  establishing  a  Legendre 
polynomial  expansion  that  approximated  the  anamorphic  wavefronts  pre¬ 
viously  designed,  as  these  were  known  to  represent  reasonable  start- 
i  ihg  points  for  combiner  design  optimization.  We  used  the  fourth- 
order  expansion  outlined  in  the  previous  section  to  establish  a 
spherical  wave  that  diverged  from  (or  converged  to)  the  vertical 
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point  focus  for  each  construction  wavefront.  Because  we  assumed  a 
curved  hologram  substrate  (with  a  radius  of  curvature  of  550  mm), 
the  fourth-order  expansion,  which  was  derived  for  flat  holograms, 
did  not  result  in  a  combiner  with  the  desired  focal  length.  There¬ 
fore,  working  with  one  wavefront  at  a  time,  we  scaled  xQ  and  zQ 
to  adjust  the  coefficients  in  Eq.  (45)  until  the  vertical  focal 
length  of  the  hologram  (191  mm)  was  the  same  as  that  of  the  conven¬ 
tional  holographic  combiner  described  in  Section  1 1 1 . 1 . 

The  resulting  combiner  had  the  desired  focal  length,  but  required 
an  intermediate  image  with  large  amounts  of  tilt  and  astigmatism. 
We  found  that  manual  adjustment  of  certain  weighting  coefficients 
was  effective  in  improving  different  aspects  of  the  hologram  perfor¬ 
mance.  For  example,  the  primary  curvature  of  each  construction  wave- 
front  is  controlled  by  the  second-order  weighting  coefficients  C?n 
and  Cq2,  which  determine  the  strength  of  the  x  and  y  phase 
variations,  respectively.  These  coefficients  were  therefore  effec¬ 
tive  at  correcting  the  hologram  astigmatism.  Similarly,  we  observed 
that  the  C3Q  coefficient  was  useful  for  adjusting  the  tilt  re¬ 
quired  of  the  image  surface  over  the  vertical  field  whereas  C4Q 
was  effective  in  adjusting  the  curvature  required  of  the  image  sur¬ 
face  over  the  vertical  field.  Because  of  the  symmetry  of  the  optical 
system  in  the  horizontal,  or  y,  direction,  the  set  of  available 
coefficients  was  limited  to  those  that  weighted  polynomials  that  were 
even  in  y,  such  as  C1Q,  C12,  C14,  C^g,  etc. 

We  found  that  by  careful  adjustment  of  the  coefficients,  we  were 
able  to  establish  construction  beam  wavefronts  that  resulted  in  good 
hologram  performance.  The  tilt  required  of  the  intermediate  image 
surface  was  reduced  and  the  hologram  astigmatism  was  corrected  over 
most  of  the  f ield-of-view.  The  distortion  introduced  by  the  hologram 
was  significantly  reduced  over  that  of  the  baseline  hologram  design 
of  Section  1 1 1 . 1 .  The  major  remaining  aberration  introduced  by  the 
hologram  was  coma,  with  a  relatively  small  amount  of  spherical 
aberration . 
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We  completed  the  aspheric  combiner  design  by  a  computer  optimiza¬ 
tion  of  the  design  that  we  had  arrived  at  through  manual  adjustment 
of  the  wavefront  weighting  coefficients.  Although  we  had  hoped  for 
considerable  reductions  in  coma,  we  found  that  the  optimization  did 
not  result  in  significant  overall  performance  improvements  of  the 
combiner . 

The  completed  aspheric  combiner  design  had  a  focal  length  of  200 
mm  with  construction  wavefronts  described  by  defining  the  wavefront 
phase  on  the  hologram  surface  in  terms  of  Legendre  polynomials.  The 
weighting  coefficients  for  both  wavefronts  are  given  in  Table  1. 
The  very  small  values  of  several  of  the  reference  beam  coefficients 
are  not  significant.  They  resulted  from  the  fourth-order  expansion 
used  to  establish  the  initial  construction  beam  parameters,  and  were 
not  changed  thereafter.  Figures  8  and  9  show  ray  traces  of  the  ref¬ 
erence  and  object  wavefronts,  respectively;  a  considerable  amount  of 
anamorphism  is  apparent  in  both  wavefronts. 

C.  ASPHERIC  COMBINER  PERFORMANCE 

In  the  HUD  optical  system,  the  relay  lens  forms  an  intermediate, 
real  image  of  the  CRT  display  at  the  focal  surface  of  the  holographic 
combiner,  so  that  the  latter,  in  turn,  will  collimate  the  display 
illumination.  Because  of  the  aberrations  and  distortions  introduced 
by  the  hologram,  the  relay  lens  must  form  an  intermediate  image 
aberrated  and  distorted  in  such  a  way  as  to  exactly  compensate  for 
the  degradation  caused  by  the  hologram.  We  can  therefore  gauge  the 
hologram  performance  by  examining  the  intermediate  image  required  by 
the  hologram  in  order  to  form  a  high  quality  image  when  viewed  from 
the  exit  pupil.  A  convenient  way  to  ascertain  the  intermediate  image 
characteristics  is  to  trace  collimated  ray  bundles  at  several  field 
angles  from  the  exit  pupil  to  the  hologram  and  determine  the  image 
formed  by  the  hologram.  The  image  so  formed  represents  exactly  the 
image  that  the  relay  lens  will  be  required  to  form. 
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TABLE  1 

WEIGHTING  COEFFICIENTS  FOR  ASPHERIC  COMBINER 
CONSTRUCTION  WAVEFRONTS 


Object  Beam 


C10 

-258580.0 

C20 

-17417.0 

C30 

6000.0 

o 

O 

-2500.0 

C02 

-39000.0 

C12 

5444.0 

C22 

511.6 

O 

o 

3000.0 

C24 

3000.0 

C34 

6000.0 

Reference  Beam 


C10 

-263890.0 

C20 

1750.0 

o 

on 

o 

17.0 

C40 

0.123 

C02 

3248.0 

C12 

55.27 

C22 

-0.523 

C04 

-0.175 
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(b)  Horizontal  Ray  Trace 


Figure  8.  Ray  Traces  of  the  Aspheric  Combiner  Reference 
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(a)  Vertical  Ray  Trace 


Figure  9.  Ray  Traces  of  the  Aspheric  Combiner  Object  Beam 
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To  determine  the  shape  and  tilt  of  the  intermediate  image  sur¬ 
face,  we  traced  a  small  pencil  of  parallel  rays  at  a  number  of  field 
angles  from  the  center  of  the  exit  pupil  to  the  intermediate  image. 
Figure  10  shows  plots  of  the  image  locations  formed  by  the  ray  pen¬ 
cils  for  both  elevation  and  azimuth  variations  about  the  center  of 
the  field.  In  both  plots,  the  horizontal  axis  is  coincident  with 
the  optical  axis  of  the  system  (which  coincides  with  the  chief  ray 
in  the  center  of  the  f ield-of-view) .  In  the  plots,  the  hologram  lies 
to  the  left  of  the  image  surface,  and  the  labels  on  the  horizontal 
axis  represent  distance  from  the  center  of  the  hologram.  In  Fig¬ 
ure  10a,  the  image  surface  is  tilted  by  41°  from  the  vertical,  and 
it  is  slightly  "s"-shaped.  The  image  surface  in  azimuth  is  slightly 
curved,  and  is  convex  toward  the  hologram.  The  images  formed  by 
vertical  and  horizontal  fans  are  indicated  by  the  two  symbol  types 
on  the  plots;  the  astigmatism  correction  is  good  at  most  field 
angles. 

We  plotted  ray  intercept  curves  as  shown  in  Figure  11  to  deter¬ 
mine  the  amount  and  types  of  aberrations  introduced  by  the  combiner. 
Data  for  a  ray  intercept  plot  are  obtained  by  tracing  a  fan  of  paral¬ 
lel  rays  from  the  exit  pupil  to  the  image  formed  by  the  hologram, 
and  at  the  image  plane  calculating  the  x  and  y  (vertical  and  horizon¬ 
tal)  displacements  between  each  ray  and  the  chief  ray  (the  chief  ray 
is  the  ray  that  passes  through  the  center  of  the  exit  pupil).  The 
image  plane  intercepts  of  the  individual  rays  relative  to  the  chief 
ray  intercept  are  plotted  as  a  function  of  the  intercepts  of  the  rays 
in  the  input  fan  (at  the  exit  pupil  of  the  HUD  system).  For  a  per¬ 
fect  image,  all  points  would  lie  on  the  horizontal  axis  as  the  image 
plane  intercepts  would  all  be  zero  relative  to  the  chief  ray 
intercept . 

Figure  11  shows  five  sets  of  ray  intercept  curves  made  for  five 
different  field  angles.  At  each  field  angle,  two  plots  are  drawn 
corresponding  to  vertical  and  horizontal  fans  passing  through  the 
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(a)  Image  Surfaces  in  the  Vertical  Plane  Formed  by  Variation 
of  the  Vertical  Field  Angle  Over  a  12°  Field.  The  Symbols 
Correspond  to  an  Angular  Increment  of  2.0° 


(b)  Image  Surfaces  in  the  Horizontal  Plane  Formed  by  Voriotion 
of  the  Horizontal  Field  Angle  Over  a  30°  Field.  The  Symbols 
Correspond  to  an  Angular  Increment  of  2.5° 


Figure  10.  Image  Surfaces  Formed  by  the  Aspheric  Combiner 
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Figure  11.  Ray  Intercept  Curves  for  the  Aspheric  Combiner  at  Five  Vertical  Field  Angles. 

The  distance  between  tick  marks  on  the  horizontal  axes  is  10  mm;  the  distance 
between  tick  marks  on  the  vertical  axes  is  1.0  mm. 


exit  pupil.  The  tick  marks  along  the  vertical  axes  each  represent 
an  interval  of  1.0  mm  at  the  image  plane  whereas  the  tick  marks  along 
the  horizontal  axes  represent  intervals  of  10  mm  in  the  input  fan. 
The  vertical  and  horizontal  input  fans  usually  extend  across  the 
vertical  and  horizontal  diameters  of  the  exit  pupil,  respectively. 
The  curves  correspond  to  elevation  angles  (denoted  by  a)  of  0,  -3, 
-6,  -9,  and  -12°;  the  azimuth  angle  in  each  case  is  0°.  The  solid 
curves  in  Figure  11  represent  the  vertical  intercepts  of  the  rays 
and  the  dashed  curves  represent  horizontal  intercepts.  For  a  ho^- 
gram  with  a  200  mm  focal  length,  an  intercept  of  1  mm  corresponds  to 
5  mr  of  ray  error  at  the  exit  pupil. 

In  Figure  11,  the  solid  curves  are  generally  concave  downward, 
indicating  coma  in  the  image.  The  dashed  curves  tend  to  be  slightly 
"s"-shaped,  which  indicates  the  presence  of  spherical  aberration. 
Ray  errors  at  the  edges  of  the  exit  pupil  are  on  the  order  of  10  to 
15  mr. 

To  obtain  a  qualitative  measure  of  distortion,  we  traced  a  set 
of  rays  from  the  center  of  the  exit  pupil  to  a  plane  tangent  to  the 
center  of  the  intermediate  image  surface.  The  rayset  comprised  rays 
distributed  uniformly  over  the  instantaneous  f ield-of-view.  Figure 
12  shows  the  intercepts  of  the  rayset  at  the  intermediate  image 
plane;  keystone  distortion  is  evident  with  points  of  constant  eleva¬ 
tion  lying  on  lines  curved  upward. 
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Figure  12.  Ray  Intercepts  Formed  by  the  Aspheric  Combiner  at  the  Intermediate  Image  Plane 
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SECTION  IV 
SYSTEM  DESIGN 


The  system  design  involved  the  addition  of  the  fold  mirror,  relay 
lens,  and  CRT  image  plane  to  the  aspheric  combiner  design  established 
previously.  Most  of  the  design  effort  concentrated  on  computer  op¬ 
timization  of  the  system  parameters,  with  the  goal  of  minimizing 
aberrations  and  distortion  while  maintaining  high  diffraction  effi¬ 
ciency  at  the  combiner. 

1.  SYSTEM  COMPONENTS 

System  components  added  were  the  fold  mirror,  relay  lens,  CRT 
image  plane,  and  ultimately,  the  CRT  faceplate.  The  fold  mirror 
folds  the  optical  path  between  the  F-16  optical  module  and  the  com¬ 
biner.  The  mirror  was  made  to  be  slightly  convex  to  reduce  the  angu¬ 
lar  spread  of  the  ray  bundles  between  the  relay  lens  and  the  mirror. 

The  initial  relay  lens  was  in  triplet  form,  and  comprised  a  doub¬ 
let  of  positive  lens  elements  (nearest  the  fold  mirror),  a  single 
negative  element,  and  a  triplet  of  positive  elements.  The  six  ele¬ 
ments  were  initially  coaxial,  and  provided  a  magnification  close  to 
uni ty. 

The  CRT  image  plane  was  initially  located  about  85  mm  from  the 
nearest  surface  of  the  relay  lens,  and  its  vertical  extent  corre¬ 
sponding  to  the  instantaneous  vertical  field  was  37  mm. 

2.  DESIGN  PROCEDURE 

Because  of  the  large  number  of  parameters  describing  the  system, 
the  design  procedure  relied  on  simultaneous  multivariable  optimiza¬ 
tion  with  the  HOAD  program.  The  optimization  algorithm  was  based  on 
the  damped  least  squares  method  that  is  commonly  used  in  optical  de¬ 
sign  as  it  is  well  suited  to  optimization  of  nonlinear  systems.  The 
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optimization  process  reduces  a  merit  function  that  is  the  sum  of  the 
squat  s  of  individual  system  errors  or  defects. 


In  our  design  optimization,  the  merit  function  comprised  three 
groups  of  system  defects:  r'ay  errors,  distortion  errors,  and  effi¬ 
ciency  errors.  Ray  errors,  or  aberrations,  were  measured  at  several 
exit  pupil  locations  and  at  several  field  angles.  For  each  field 
angle,  we  traced  a  set  of  5  or  6  parallel  rays  from  the  exit  pupil 
to  the  CRT  image  plane;  the  rayset  included  a  chief  ray  passing 
through  the  center  of  the  exit  pupil  with  the  remaining  rays  distrib¬ 
uted  along  the  vertical  and  horizontal  centerlines  of  the  exit  pupil. 
The  ray  errors  were  calculated  as  the  image  plane  intercepts  of  the 
individual  rays  relative  to  the  chief  ray  intercept. 

Distortion  errors  were  measured  by  tracing  a  set  of  rays,  dis¬ 
tributed  in  field  angle  over  the  instantaneous  f ield-of-view,  from 
the  center  of  the  exit  pupil  to  the  CRT  image  plane.  In  a 
distortion-free  system,  the  rayset  would  form  a  rectangular  grid  at 
the  CRT.  The  distortion  contribution  to  the  merit  function  was  made 
up  of  departures  of  the  actual  ray  intercepts  from  a  perfectly  rec¬ 
tangular  grid,  where  the  grid  spacing  was  determined  as  the  distance 
between  the  rays  nearest  to  the  central  ray.  To  simplify  the  distor¬ 
tion  measurements,  we  compared  distances  at  the  CRT  that  corresponded 
to  equal  increments  in  field  angle.  Consequently,  the  optimization 
was  based  on  satisfying  the  equation  h  =  ko  where  h  is  the  image 
height  at  the  CRT,  9  is  the  field  angle  at  the  exit  pupil,  and  k  is 
a  constant;  both  h  and  e  are  measured  relative  to  the  center  of  the 
f ield-of-view.  Since  we  anticipated  that  the  ultimate  display  source 
would  be  a  matrix  addressed  device,  for  most  of  our  design  investiga¬ 
tions,  we  allowed  the  constant  to  be  different  in  the  horizontal  and 
vertical  directions. 

Five  rays  were  used  to  monitor  the  hologram  efficiency,  primarily 
to  prevent  the  design  parameters  from  drifting  away  from  the  geometry 
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required  to  maintain  high  diffraction  efficiency  over  both  the  extent 
of  the  exit  pupil  and  the  f ield-cf-view.  The  coupled  wave  analysis 
[6]  used  to  calculate  the  hologram  diffraction  efficiency  is  correct 
only  for  rays  lying  in  a  plane  of  incidence  that  also  contains  the 
grating  vector  (a  vector  that  is  perpendicular  to  the  grating  planes 
within  the  hologram  medium);  in  the  HUD  system,  this  is  the  vertical 
plane  lying  in  the  center  of  aircraft  (the  xz-plane  in  our  coordinate 
frame).  Consequently,  the  five  rays  used  to  sample  the  hologram 
efficiency  all  lay  in  the  xz-plane.  Three  rays  passed  through  the 
center  of  the  exit  pupil  at  the  top,  center,  and  bottom  of  the  IFOV 
(0°,  -6°,  and  -12°,  respectively).  Of  the  two  additional  rays,  one 
passed  through  the  top  of  the  exit  pupil  at  an  angle  corresponding 
to  the  bottom  of  the  TFOV  (-16°)  and  the  other  passed  through  the 
bottom  of  the  exit  pupil  at  the  top  of  the  TFOV  (+4°).  For  the  effi¬ 
ciency  calculations,  we  assumed  a  hologram  thickness  of  10  um  and  a 
refractive  index  modulation  of  0.04. 

We  changed  the  rayset  occasionally  as  the  design  evolved.  The 
initial  optimization  was  performed  with  rays  covering  an  undersized 
exit  pupil  and  with  coarse  sampling  for  distortion  correction.  As 
the  design  improved,  the  rayset  was  redistributed  to  correspond  to 
the  full  exit  pupil  size  and  rays  were  added  for  increased  distortion 
sampling.  The  rayset  used  in  the  final  optimization  had  a  total  of 
77  rays,  of  which  37  were  used  to  determine  ray  errors,  45  were  used 
i.o  measure  distortion,  and  5  were  used  to  measure  the  hologram  effi¬ 
ciency  l some  rays  contributed  to  more  than  one  error  calculation). 
In  addition  to  changing  the  rayset,  we  also  altered  as  appropriate 
the  weighting  with  which  the  squares  of  the  individual  errors  were 
added  to  the  merit  function. 

The  HUD  optical  system  was  described  by  more  than  100  parameters 
that  were  reasonable  optimization  variables.  Of  these  parameters, 
64  were  construction  beam  weighting  coefficients  while  most  of  the 
remaining  variables  were  associated  with  the  relay  lens.  Although 


we  did  not  optimize  with  all  of  the  variables  at  once,  we  did  use 
most  of  the  variables  over  the  course  of  the  investigation,  with  a 
typical  optimization  run  including  25  variables  selected  from  the 
total  set.  We  found  that  occasionally  adding  a  new  variable  (such 
as  a  lens  decentration  or  tilt)  was  effective  at  increasing  the  rate 
of  improvement  in  the  optimization  results. 

3.  DESIGN  RESULTS 

Computer  optimization  was  carried  out  using  the  HOAD  program. 
The  optimization  resulted  in  an  aspheric  combiner  that  differed  con¬ 
siderably  from  the  initial  design  and  in  a  relay  lens  with  tilted 
and  decentered  elements. 

a.  ASPHERIC  COMBINER 

The  resultant  aspheric  combiner  was  defined  by  a  greater  number 
of  construction  beam  weighting  parameters  than  had  been  used  in  the 
initial  combiner  design.  As  in  the  initial  design,  both  construction 
wavefronts  were  anamorphic ;  the  coma,  however,  which  had  been  ob¬ 
served  in  the  object  beam  moved  to  the  reference  beam.  Ray  traces 
of  the  reference  and  object  beams  are  shown  in  Figures  13  and  14, 
respectively.  In  addition,  the  radius  of  curvature  of  the  combiner 
decreased  slightly  from  550.1  mm  to  549,3  mm.  The  weighting  coeffi¬ 
cients  that  define  the  reference  and  object  beams  are  shown  in 
Table  2.  Whereas  a  total  of  18  coefficients  had  been  used  to  specify 
the  construction  beams  of  the  initial  combiner  design,  the  construc¬ 
tion  beams  of  the  final  combiner  design  utilized  a  total  of  66 
coefficients. 

The  performance  of  the  aspheric  combiner  is  indicated  by  the  ray 
intercept  curves  in  Figure  15.  The  coma  that  was  evident  in  the 
initial  combiner  design  (Figure  11)  was  decreased  but  not  eliminated, 
and  some  astigmatism  has  been  added.  The  astigmatism  compensates 
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(a)  Vertical  Ray  Trace 
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Figure  13.  Ray  Traces  o<  the 


Final  Aspheric  Combiner  Reference  Beam 
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Figure  14.  Ray  Traces  of  the  Final  Aspheric  Combiner  Object  Beam 
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TABLE  2 

WEIGHTING  COEFFICIENTS  FOR  THE  FINAL  ASPHERIC 
COMBINER  CONSTRUCTION  WAVEFRONTS 


Object  Beam 


CIO 

-258580.0 

C04 

150.60 

C20 

-9133.5 

C14 

8776.4 

C30 

4621  .4 

C24 

1622.9 

C40 

-1627.9 

-34 

2711.3 

C50 

-301.17 

C44 

122.31 

C60 

-239.70 

C70 

620.91 

CO  6 

358.37 

C80 

230.79 

C16 

2383.9 

C90 

526.96 

C26 

-23.742 

C02 

-35240.0 

C08 

62.332 

C12 

31463.0 

C18 

57.548 

C22 

3310.2 

C32 

3618.2 

C42 

153.54 

C52 

-2303.4 

C62 

-177.46 

Reference  Beam 

CIO 

-262889.0 

C04 

-1730.6 

C20 

-910.09 

Cl  4 

-1578.0 

C30 

-2515.5 

C24 

-180.31 

C40 

752.01 

C34 

-507.19 

C50 

-1033.5 

C44 

146.00 

C60 

-417.73 

C54 

-283.96 

C70 

-683.85 

C64 

-370.24 

C80 

-85.929 

C74 

354.98 

C90 

-316.70 

C84 

-494.25 

C02 

-43.005 

C06 

249.29 

C12 

2863.4 

C16 

-179.36 

C22 

1346.7 

C26 

55.215 

C32 

6628.0 

C36 

375.31 

C42 

-269.56 

C46 

-365.22 

C52 

-68.542 

C56 

202.08 

C62 

-165.57 

C66 

-542.16 

C72 

-141.01 

C82 

-156.42 

C08 

-20.738 

C92 

-123.20 

C18 

-299.30 

C28 

-193.14 

C38 

-174.55 

C48 

10.498 

49 


for  astigmatism  introduced  by  the  tilts  and  decenters  of  the  relay 
lens  that  were  required  to  correct  for  the  coma. 

b.  RELAY  LENS 

The  relay  lens  retained  its  basic  triplet  form  except  that  during 
the  design  optimization,  we  split  the  center  negative  element  into 
two  negative  elements.  In  the  course  of  the  optimization  procedure, 
one  of  the  negative  elements  became  a  positive  element,  with  the 
second  element  taking  on  stronger  negative  power. 

The  design  of  the  relay  lens  proceeded  with  consideration  given 
to  lens  manufacture  and  assembly.  All  surfaces  were  kept  spherical, 
and  the  number  of  tilts  was  controlled  by  arranging  the  lens  elements 
into  three  groups. 

Figure  16  shows  the  configuration  of  the  relay  lens  that  resulted 
from  the  system  optimization.  In  Figure  16,  the  broken  line  connects 
the  vertices  of  the  lens  elements  and  the  solid  line  represents  a 
ray  traced  from  the  center  of  the  exit  pupil  at  the  central  field 
angle.  The  lens  elements  are  arranged  in  three  groups  that  are 
tilted  and  decentered  with  respect  to  one  another.  The  group  nearest 
the  fold  mirror  is  a  positive  doublet  in  which  both  elements  are 
coaxial.  The  second,  or  center  group,  comprises  a  positive  and  a 
negative  element  which  are  decentered  vertically  with  respect  to  each 
other,  but  which  have  parallel  axes  of  rotation.  The  third  group  is 
a  positive  triplet  in  which  the  element  nearest  to  the  CRT  is  decen¬ 
tered  vertically  with  respect)  to  the  first  two  elements;  its  axis  is 
parallel  to  the  axis  of  the  fjirst  two  elements. 

The  relay  lens  prescription  is  given  by  group  in  Table  3,  in 
which  a  positive  radius  of  curvature  indicates  that  the  center  of 
curvature  lies  to  the  left  of  ithe  surface  as  shown  in  Figure  16,  and 
a  positive  decentration  indicates  a  displacement  upward  in  Figure  16. 
The  figure  also  indicates  the  surface  and  group  numbering. 


5' 


TABLE  3 

RELAY  LENS  PARAMETERS  BY  GROUP 


Group  1 
Surface  No. 

Glass  Type 

Index 

Thickness 

(mm) 

Radius 

(mm) 

1 

LaSF  N3 

1.8127 

20.000 

86.461 

2 

AIR 

1.00 

1 .300 

154.417 

3 

8K  7 

1.5188 

45.500 

94.482 

4 

AIR 

1.00 

-193.701 

Group  2 

1 

BK  7 

1.5188 

7.083 

57.433 

2 

AIR 

1.00 

13.303 

227.443 

3 

8K  7 

1.5188 

5.000 

-44.272 

4 

AIR 

1.00 

63.886 

Group  3 

1 

BaF  8 

1.6269 

21.330 

-768.038 

2 

AIR 

1.00 

1.000 

-116.258 

3 

LaSF  N3 

1.8127 

14.000 

326.690 

4 

AIR 

1.00 

1.000 

-100.030 

5 

SF  1 

1.7231 

19.998 

83.182 

6 

AIR 

1.00 

-15028.439 

Decentration 


-3.341 


-8.604 
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Table  4  gives  the  separation  and  tilt  parameters  between  the  re¬ 
lay  lens  groups  and  between  the  relay  lens  and  the  CRT  faceplate. 
The  location  of  the  vertex  of  the  first  surface  of  group  2  with  re¬ 
spect  to  group  1  is  found  as  follows.  Proceed  1  mm  along  the  axis 
of  group  1  as  it  leaves  the  vertex  of  surface  4;  move  5.568  mm  per¬ 
pendicular  to  the  axis  in  an  upward  direction.  Then  rotate  8.827° 
to  establish  the  direction  of  the  axis  of  group  2,  where  a  posit  _■ 
rotation  corresponds  to  rotating  the  axis  clockwise  in  Figure  16. 

c.  SYSTEM  CONFIGURATION 

Figure  17  shows  the  system  configuration.  The  parameters  that 
locate  the  combiner,  fold  mirror,  and  first  surface  of  the  relay  lens 
are  given  in  Table  5,  in  which  the  separation  distance  is  the  dis¬ 
tance  between  the  element  centers  and  the  separation  angle  is  the 
angle  between  the  surface  normal  and  the  line  connecting  the  element 
centers.  For  example,  the  distance  between  the  center  of  the  exit 
pupil  and  combiner  is  625  mm,  measured  along  a  line  that  makes  an 
angle  of  6°  to  the  normal  that  passes  through  the  center  of  the  exit 
pupil.  The  exit  pupil  is  assumed  to  be  vertical  with  respect  to  the 
aircraft  cockpit.  The  sign  convention  for  the  tilts  in  Table  5  is 
the  same  as  for  the  group  tilts  in  Table  4.  The  fold  mirror  is  con¬ 
vex  with  a  radius  of  curvature  of  749.38  mm,  and  we  assumed  a  BK  7 
faceplate  for  the  CRT  with  a  thickness  of  3.2  mm. 

Figure  18  shows  horizontal  and  vertical  ray  traces  through  the 
system,  and  indicates  how  the  system  fits  within  the  aircraft  cock¬ 
pit.  In  Figure  18a,  the  line  labelled  "CLEARANCE  LINE"  is  a  line 
0.75  in.  inside  the  inner  surface  of  the  F— 1 6  canopy,  which  is  shown 
without  a  label.  Although  the  uppermost  ray  in  Figure  18a  crosses 
the  clearance  line,  that  ray  is  outside  the  ray  envelope  required  by 
the  field-of-view  specifications  and  is  shown  for  completeness  in 
the  ray  trace. 


TABLE  4 

RELATIVE  LOCATION  OF  RELAY  LENS  GROUPS 


Axial  Separation 
(mm) 

Decemration 

(mm) 

Group  1-2 

1.000 

5.568 

Group  2-3 

11.103 

0.200 

Group  3-CRT 

52.024 

0.000 

TABLE  5 

RELATIVE  LOCATION  OF  SYSTEM 

COMPONENTS 

Separation 

Oistance 

(mm) 

Separation 

Angle 

(deg) 

Tilt 

(deg) 

Exit  Pupil-to- 
Combiner 

625.00 

6.00 

24.00 

Combiner-to- 
Fold  Mirror 

117.07 

32.65 

-4.37 

Fold  Mirror-to- 
Relay  Lens 

136.55 

35.59 

-39.37 
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SECTION  V 

SYSTEM  PERFORMANCE 

The  system  was  character ized  primarily  in  terms  of  distortion, 
ray  errors,  and  hologram  efficiency.  The  system  met  the  goal  of  less 
than  0.5  percent  distortion  over  most  of  the  instantaneous  field-of- 
view,  and  the  hologram  efficiency  was  also  satisfactory  over  most  of 
the  field.  The  ray  errors  exceeded  the  design  goals  over  much  of 
the  f ield-of-view  and  exit  pupil,  but  were  generally  not  excessive. 

1.  DISTORTION 

Distortion  was  measured  quantitatively  by  tracing  rays  from  the 
center  of  the  exit  pupil  at  different  field  angles  and  finding  their 
intercepts  at  the  CRT  image  plane.  Figure  19  shows  plots  of  distor¬ 
tion  as  a  function  of  field  angle  over  the  instantaneous  vertical 
and  horizontal  f i eld-of-view.  For  these  plots,  the  distortion  was 
measured  as  follows.  A  measurement  was  made  very  close  to  the  center 
of  the  field  to  determine  aS/ao,  the  variation  aS  in  image  intercept 
with  a  small  variation  as  in  field  angle,  from  which  we  calculated  a 
predicted  image  intercept  Sp  (relative  to  the  intercept  at  the 
center  of  the  field)  as  a  function  of  the  field  angle  9, 


We  calculated  the  percent  distortion  D  by  comparing  the  measured  S 
values  against  the  predicted  Sp  values,  as 

S  -  S 

D  =  -E- -  x  1 00 

P 

In  the  horizontal  direction,  the  distortion  meets  the  goal  of 
0.5  percent  over  the  instantaneous  and  total  fields-of-view  of  ±10°, 
with  a  peak  distortion  of  0.46  percent  distortion  at  ±10°.  In  the 
vertical  direction,  the  distortion  meets  the  goal  from  -12°  to  nearly 
-2°,  increasing  to  a  peak  of  2.8  percent  at  an  elevation  angle  of  0° . 
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Figure  19.  Distortion  Over  the  Instantaneous  Fieid-of-View 
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As  the  distortion  was  evaluated  over  the  instantaneous  field  only 
during  the  design  optimization,  we  only  show  distortion  performance 
over  the  instantaneous  field.  This  was  a  result  of  measuring  the 
distortion  from  the  center  of  the  exit  pupil  while  defining  the  com¬ 
biner  size  such  that  field  angles  outside  the  (vertical)  instantane¬ 
ous  field  require  head  motion  away  from  the  exit  pupil  center.  Al¬ 
though  the  hologram  could  have  been  extended  beyond  the  intended 
limits  and  rays  traced  from  the  center  of  the  exit  pupil  over  the 
total  vertical  field,  such  measurements  would  be  of  limited  meaning. 
Perhaps  a  more  appropriate  solution  would  have  been  to  define  the 
distortion  relative  to  rays  traced  from  locations  in  the  exit  pupil 
that  varied  as  a  function  of  field  angle  so  that  the  total  field 
could  have  been  appropriately  sampled  and  measured. 

Figure  20  shows  the  image  plane  intercepts  of  a  set  of  rays 
traced  from  the  center  of  the  exit  pupil  and  distributed  over  the 
instantaneous  f ield-of-view  in  2°  increments.  The  grid  formed  by 
the  intercepts  has  relatively  straight  lines  and  regular  spacing, 
although  some  distortion  is  still  evident,  particularly  at  the  cor¬ 
ners  of  the  grid. 

The  distortion  performance  represents  a  significant  improvement 
over  that  of  previous  holographic  HUD  systems  that  we  have  investi¬ 
gated.  Some  aspects,  however,  of  the  distortion  performance  of  the 
conventional  holographic  systems  were  comparable  to  the  corresponding 
performance  of  the  aspheric  holographic  system,  but  these  aspects 
differed  with  the  individual  designs  and  they  represented  only  a 
portion  of  the  distortion  performance.  Figure  21  shows  curves  of 
the  distortion  magnitude  as  functions  of  horizontal  and  vertical 
field  angles  for  the  current  system,  and  for  an  average  of  data  from 
conventional  holographic  HUD  systems. 


The  distortion  improvement  is  particularly  evident  when  plots  of 
image  plane  intercepts  (as  in  Figure  20)  are  compared.  Figure  22, 
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Figure  20.  Image  Plane  Intercepts  of  a  Rayset  Distributed  Over  the  Instantaneous 
Field-of-View  in  2°  Increments 
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Conventional  system  distortion  shown 
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System  Distortion  Magnitude  Compared  with  Average  Distortion  Magnitude  of  Conventional  Holographic  HUD  Systems. 
Conventional  system  distortion  shown  dashed. 
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Figure  22.  Image  Plane  Intercepts  of  a  Rayset  Uniformly  Distributed  Over  a  12  x  20° 
Field-of-View  for  a  Conventional  Holographic  HUD  System 
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for  example,  shows  the  image  plane  intercepts  corresponding  to  a  12  x 
20°  field  for  a  system  in  which  the  preliminary  combiner  design  dis¬ 
cussed  in  Section  III. lb  was  combined  with  the  initial  relay  lens 
design.  While  the  system  was  far  from  optimum,  the  distortion  shown 
in  Figure  22  is  indicative  of  the  distortion  typically  found  in  holo¬ 
graphic  HUD  systems. 

2.  RAY  ERRORS 

We  measured  ray  errors  by  generating  ray  intercept  curves  at  the 
CRT  image  surface  and  by  measuring  the  reticle  parallax  and  reticle 
accuracy  at  the  exit  pupil.  Ray  intercept  curves  are  shown  in  Fig¬ 
ure  23  for  angles  covering  the  instantaneous  field-of-view.  These 
curves  were  generated  by  tracing  fans  of  parallel  rays  from  the  exit 
pupil  to  the  CRT  image  surface  as  described  previously  for  the  inter¬ 
mediate  image  surface.  In  Figure  23,  the  vertical  axes  correspond 
to  a  displacement  of  ±1  mm,  which,  in  turn,  corresponds  to  an  angular 
error  at  the  exit  pupil  of  approximately  ±6.8  mr.  The  curves  repre¬ 
sent  elevation  angles  (denoted  by  a)  of  0,  -3,  -6,  -9,  and  -12°,  and 
azimuth  angles  (denoted  by  b)  of  0,  -3.3,  -6.7,  and  -10°;  because  of 
symmetry,  the  curves  would  have  the  same  appearance  for  positive 
azimuth  angles. 

The  ray  intercept  curves  are  useful  as  they  indicate  the  types 
of  aberration  remaining  in  the  optical  system.  At  the  field  angles 
shown,  sphrical  aberration  is  the  primary  residual  aberration,  with 
a  lesser  amount  of  uncorrected  coma  present. 

Figure  24  shows  ray  intercept  curves  made  at  elevation  angles  of 
+3°  and  -15°,  3°  in  each  direction  beyond  the  limits  of  the  vertical 
instantaneous  field-of-view.  For  the  curves  at  +3°,  the  fans  were 
centered  1.0  cm  below  the  center  of  the  eyebox  whereas  for  the  curves 
at  -15°,  the  fans  were  centered  1.0  cm  above  the  center  of  the  eye- 
box.  Considerable  coma  and  spherical  aberration  are  evident  in  the 
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a  *  0.00 


Vertical  Fan 


Horizontal  Fan 


(b)  Horizontal  Field  Angles 


Figure  23.  Ray  Intercept  Curves  for  the  Optical  System  at  Five  Vertical  Field  Angles  and  Three 
Horizontal  Field  Angles.  The  distance  between  tick  marks  on  the  horizontal  axes  is 
10  nun;  the  distance  between  tick  marks  on  the  vertical  axes  is  0.5  mm. 
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Figure  24.  Ray  Intercept  Curves  for  the  Optical  System  at  Two  Vertical  Field  Angles  Beyond 
the  Extent  of  the  Vertical  Instantaneous  Field-of-View.  The  distance  between 
tick  marks  on  the  horizontal  axes  is  10  mm;  the  distance  between  tick  marks  on 
the  vertical  axes  is  0.5  mm. 
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curves  drawn  for  the  horizontal  fans,  a  result  of  our  limiting  the 
rayset  used  in  the  system  optimization  to  the  instantaneous  field. 

Although  the  ray  intercept  curves  provide  useful  insights  into 
the  behavior  and  performance  of  the  system,  the  performance  criteria 
of  greatest  ultimate  interest  are  the  accuracy  and  parallax  errors 
at  the  exit  pupil.  As  a  measure  of  these  errors,  we  determined  the 
reticle  accuracy  and  the  reticle  parallax  as  follows.  In  each  case, 
the  method  was  based  on  finding  the  actual  ray  angles  at  the  two  eye 
positions  that  corresponded  to  some  nominal  field  angle.  Assuming  a 
distortion-free  system,  we  determined  a  CRT  location  corresponding 
to  a  nominal  field  angle.  Then,  for  a  given  head  location,  we  found 
the  rays  passing  through  the  two  eye  locations  that  reached  the  CRT 
location  corresponding  to  the  nominal  field  angle.  The  ray  errors 
were  the  differences  in  elevation  a  a  and  in  azimuth  a  b  between  the 
actual  ray  angles  and  the  nominal  field  angles.  The  right  eye  errors 
are  thus 


ActR  =  “right  ~  “nom. 


abR  =  8right  "  Bnom. 

and  a  similar  set  of  equations  define  a«l  and  abl  for  the  left 
eye  errors. 


Reticle  accuracy  was  then  calculated  as  the  root-sum-square  of 
the  ray  errors  averaged  between  the  two  eyes,  which  can  be  expressed 
as 


Reticle  accuracy  = 


At  field  angles  and  head  locations  where  the  display  was  not  visible 
to  both  eyes,  reticle  accuracy  was  computed  as  the  root-sum-square 
of  the  errors  corresponding  to  the  eye  for  whico  the  display  was 
visible. 
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V 


Reticle  parallax  was  calculated  as  the  root-sum-square  of  the 
differences  in  the  errors  at  the  two  eye  locations,  or 


Reticle  parallax  =  [ (aor  -  &aL)^  +  ( 


4Br  -  abl 


Where  the  display  was  not  visible  to  both  eyes,  there  could  be  no 
parallax  and  therefore  no  calculation  was  made. 

We  determined  the  reticle  accuracy  and  reticle  parallax  as  a 
function  of  head  position  at  each  of  a  number  of  field  angles.  The 
size  of  the  exit  pupil,  3.0  in.  vertical  by  4.5  in.  horizontal  per¬ 
mitted  head  motion  of  ±1.5  in.  vertical  by  ±1.0  in.  horizontal,  keep¬ 
ing  both  eyes  within  the  exit  pupil.  Figure  25  shows  curves  of 
reticle  accuracy  and  reticle  parallax  as  functions  of  horizontal  head 
position  (in  inches)  at  the  center  of  the  field.  The  three  curves 
in  each  plot  represent  different  vertical  head  positions,  with  the 
solid  line  for  the  head  centered,  the  short  dashes  for  the  head 
raised  by  1.0  inch,  and  the  long  dashes  for  the  head  lowered  by  1.0 
inch.  The  data  were  taken  at  horizontal  intervals  of  0.2  inch.  The 
accuracy  data  average  to  0.88  mr  with  a  peak  of  2.1  mr  occurring  for 
the  head  raised  position.  The  parallax  errors  are  somewhat  more 
severe,  with  a  2.0  mr  average  and  a  peak  of  3.8  mr,  again  for  the 
head  raised  position.  We  found  that  over  most  of  the  f ield-of-view, 
the  parallax  errors  were  comparable  to,  or  more  severe  than,  the 
accuracy  errors. 

Similar  plots  of  reticle  accuracy  and  reticle  parallax  for  other 
field  angles  are  shown  in  Figures  26-28,  where  the  data  in  Fig¬ 
ures  26,  27,  and  28  represent  angular  displacements  of  2.5,  5.0,  and 
7.5°,  respectively,  from  the  center  of  the  f ield-of-view.  The  errors 
are  greatest  at  the  edges  of  the  exit  pupil  and  for  elevation  angles 
corresponding  to  upward  looks.  The  large  errors  at  upward  looks  are 
in  part  a  result  of  the  greater  distortion  at  those  field  angles,  as 
we  compute  the  errors  assuming  a  distortion-free  system.  Table  6 
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Reticle  Parallax  Reticle  Accuracy 

0.00  2.00  4.00  6.00  8.00  ,  0.00  2.00  4.00 


Elevation  -6.00  Azimuth  0.0 


120  -0.80  -0.40  0.00  0.40  0.80  1.20 

Horizontal  Head  Position 


Horizontal  Head  Position 


Figure  25.  Reticle  Accuracy  and  Reticle  Parallax  As  Functions  of  Head  Position  for  the 
Center  of  the  Field -of- View.  The  solid  line  corresponds  to  the  head  centered 
vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and  the  long  dashes  to 
the  head  lowered  1.0  inch. 
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Parallax  Reticle  Accuracy 


(a)  2.5°  Upward  Displacement 


Figure  26.  Reticle  Accuracy  and  Reticle  Parallax  As  Fucntions  of  Head  Position  for  Field  Angles 
Displaced  by  2.5°  from  the  Center  of  the  Field.  The  solid  line  corresponds  to  the  head 
centered  vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and  the  long  dashes  for 
the  head  lowered  1.0  inch 
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(b)  2.5°  lateral  Displacement 


Figure  26.  Reticle  Accuracy  and  Reticle  Parallax  As  Functions  of  Head  Position  for  Field  Angles 
Displaced  by  2.5°  from  the  Center  of  the  Field.  The  solid  line  corresponds  to  the  head 
centered  vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and  the  long  dashes  for 
the  head  lowered  1.0  inch. 
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Reticle  Accurocy 


Elevation  -1.00  Azimuth  0.0 


(a)  5.0°  Upward  Displacement 


Figure  27.  Reticle  Accuracy  and  Reticle  Parallax  As  Functions  of  Head  Position  for  Field 

Angles  Displaced  by  5.0°  from  the  Center  of  the  Field.  The  solid  line  corresponds 
to  the  head  centered  vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and 
the  long  dashes  for  the  head  lowered  1.0  inch. 
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Reticle  Parallax  Reticle  Accuracy 

6.00  9.00  U.00  ,0.00  2.00  4.00 
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(b)  5.0°  Lateral  Displacement 


Figure  27.  Reticle  Accuracy  and  Reticle  Parallax  As  Functions  of  Head  Position  for  Field 

Angles  Displaced  by  5.0°  from  the  Center  of  the  Field.  The  solid  line  corresponds 
to  the  head  centered  vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and 
the  long  dashes  for  the  head  lowered  1.0  inch. 
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Figure  27.  Retide  Accuracy  and  Retide  Parallax  As  Functions  of  Head  Position  for  Field 

Angles  Displaced  by  5.0“  from  the  Center  of  the  Field.  The  solid  line  corresponds 
to  the  head  centered  vertically,  the  short  dashes  to  the  head  raised  1.0  indi,  and 
the  long  dashes  for  the  head  lowered  1.0  inch 
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Reticle  Parallax  Reticle  *'rUracv 

30  12.00  18.00  24.00  ,0.00  6.00  12.00  16.00  24.00 


Horizontal  Head  Position 
(a)  7.5°  Upward  Displacement 

Figure  28.  Reticle  Accuracy  and  Reticle  Parallax  As  Functions  of  Head  Position  for  Field 

Angles  Displaced  by  7.5°  from  the  Center  of  the  Field.  The  solid  line  corresponds 
to  the  head  centered  vertically,  and  the  long  dashes  for  the  head  lowered  1.0  inch. 


81 


Reticle  Parallax  Reticle  Accuracy 

6.00  9.00  12.00  .  0.00  2.00  4.00 
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(b)  7.5°  Lateral  Displacement 

Figure  28.  Reticle  Accuracy  and  Reticle  Parallax  A s  Functions  of  Head  Postion  for  Field 

Angles  Displaced  by  7.5°  from  the  Center  of  the  Field.  The  solid  line  corresponds 
to  the  head  centered  vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and 
the  long  dashes  for  the  head  lowered  1.0  inch. 


Reticle 


(c)  7.5°  Downward  Displacement 


Figure  28.  Reticle  Accuracy  and  Reticle  Parallax  As  Functions  of  Head  Position  for  Field 

Angles  Displaced  by  7.5°  from  the  Center  of  the  Field.  The  solid  line  corresponds 
to  the  head  centered  vertically,  the  short  dashes  to  the  head  raised  1.0  inch,  and 
the  long  dashes  for  the  head  lowered  1.0  inch. 
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TABLE  6 

RETICLE  PARALLAX  AND  ACCURACY  ERRORS 


Parallax 

Accuracy 

a 

6 

Goal 

Average 

Maximum 

Goal 

Average 

Maximum 

-6.0 

0.0 

0.6 

2.01 

3.77 

0.8 

0.88 

2.10 

-3.5 

0.0 

0.6 

3.01 

8.08 

1.0 

1.39 

4.77 

-6.0 

2.5 

0.6 

3.04 

15.14 

1.0 

1.40 

8.22 

-8.5 

0.0 

0.6 

1.63 

4.51 

1.0 

0.98 

2.15 

-1.0 

0.0 

2.0 

5.04 

11.84 

2.0 

2.56 

6.06 

-6.0 

5.0 

2.0 

3.66 

10.55 

2.0 

2.04 

6.96 

■11.0 

0.0 

2.0 

2.43 

8.29 

2.0 

1.20 

3.97 

1.5 

0.0 

2.0 

7.90 

33.39 

3.0 

16.51 

35.25 

-6.0 

7.5 

2.0 

3.38 

8.26 

3.0 

2.72 

6.47 

■13.5 

0.0 

2.0 

4.94 

12.19 

3.0 

2.56 

8.17 

All  angles  in  degrees  and  errors  in  mi  1 1 iradians 
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contains  for  each  field  angle  the  maximum  and  average  error  of  the 
data  shown  in  Figures  25-28,  and  the  parallax  and  accuracy  design 
goals.  As  in  the  ray  intercept  plots,  the  elevation  angles  are  de¬ 
noted  by  a  and  the  azimuth  angles  by  6. 

3.  EFFICIENCY 

A  major  factor  affecting  the  brightness  uniformity  of  the  system 
is  the  hologram  diffraction  efficiency.  We  computed  the  diffraction 
efficiency  of  rays  distributed  over  the  vertical  field  and  over  the 
vertical  extent  of  the  exit  pupil.  All  rays  used  to  evaluate  the 
hologram  efficiency  were  in  the  vertical  plane  bisecting  the  optical 
system  (the  xz-plane  in  our  coordinate  frame),  as  the  analysis  we 
use  assumes  that  the  plane  of  incidence  is  parallel  to  the  hologram 
grating  vector. 

Figure  29  shows  curves  of  diffraction  efficiency  as  a  function 
of  vertical  field  angle  for  head  positions  corresponding  to  the  top, 
center,  and  bottom  of  the  exit  pupil.  We  assumed  a  volume  phase 
hologram  with  a  thickness  of  10  ^m  and  a  refractive  index  modulation 
of  0.04.  The  rays  sampled  the  system  at  1°  increments,  with  some 
rays  vignetted  by  the  optical  system.  From  the  center  and  bottom  of 
the  exit  pupil,  the  efficiency  exceeded  90  percent  for  all  rays. 
From  the  top  of  the  exit  pupil,  the  uniformity  is  much  less,  with 
rays  at  the  upper  field  angles  strongly  attenuated  by  the  hologram. 
If  we  assume  that  75  percent  of  the  vertical  instantaneous  field-of- 
view  should  be  visible  from  the  upper  edge  of  the  exit  pupil,  then 
the  required  field  coverage  extends  from  an  elevation  angle  of  -16° 
to  -8°;  over  this  field,  the  efficiency  is  everywhere  above  60 
percent. 

The  design  procedure  was  intended  to  produce  a  symmetric  uniform¬ 
ity,  such  that  if  the  efficiency  from  the  top  of  the  exit  pupil  de¬ 
creased  at  the  upper  elevation  angles,  then  the  efficiency  from  the 
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Figure  29.  Hologram  Efficiency  As  a  Function  of  Vertical  Field  Angle  for  Rays  Traced  from 
Three  Exit  Pupil  Locations.  The  solid  line  is  for  rays  traced  from  the  center  of  the 
pupil,  the  short  dashes  for  rays  from  the  top  of  the  pupil,  and  the  long  dashes  for 
rays  from  the  bottom  of  the  pupil. 
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bottom  of  the  exit  pupil  would  decrease  at  lower  elevation  angles. 
This  compromise  was  not  achieved  as  a  result  of  using  only  a  single 
ray  from  each  end  of  the  exit  pupil  to  sample  the  hologram  efficiency 
during  the  system  optimization.  It  would  have  been  helpful  to  have 
traced  two  rays  from  each  end  of  the  exit  pupil.  For  a  required 
vertical  instantaneous  field  of  8°  from  the  top  and  bottom  of  the 
exit  pupil,  for  example,  we  would  sample  the  hologram  efficiency  with 
rays  at  elevation  angles  of  -8°  and  -16°  from  the  top  of  the  pupil, 
and  at  angles  of  +4°  and  -4°  from  the  bottom  of  the  pupil. 

The  region  of  low  efficiency  is  limited  to  a  relatively  small 
portion  of  the  exit  pupil  and  f ield-of-view.  Figure  30  shows  effi¬ 
ciency  curves  for  exit  pupil  locations  0.75  inch  above  and  below  the 
center  of  the  pupil;  these  locations  are  half  the  distance  to  the 
limits  of  the  pupil.  The  ray  efficiencies  plotted  in  Figure  30  ex¬ 
ceed  88  percent,  and  over  most  of  the  field,  are  above  90  percent. 
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I  Diffraction  Efficiency 
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Figure  30.  Hologram  F.ffitiency  As  a  Function  of  Vertical  Field  Angle  for  Rays  Traced  from 
Two  Exit  Pupil  Locations.  The  solid  line  is  for  rays  traoed  from  0.75  inch  above 
the  pupil  center;  the  dashed  line  for  rays  traced  from  0.75  inch  below  the  pupil 
center. 


83 


SECTION  VI 

REALIZATION  OF  CONSTRUCTION  BEAMS 


The  discussion  thus  far  has  considered  the  construction  beams  to 
be  defined  simply  by  a  description  of  their  respective  phase  varia¬ 
tions  at  the  hologram.  Aside  from  the  analysis  in  Section  II  proving 
the  realizabilty  of  the  phase  functions,  we  have  not  considered  how 
the  construction  beams  would  be  implemented  to  fabricate  the  aspheric 
combiner.  In  this  section,  we  discuss  practical  considerations  for 
the  realization  of  the  analytically  defined  construction  beams. 

1.  DESIGN  PROCEDURE 


The  form  in  which  the  construction  wavefronts  are  described  does 
not  suggest  a  direct  method  for  their  generation.  A  conceptually 
straightforward  method  would  be  to  record  the  combiner  itself  as  a 
computer  generated  hologram  (CGH)  [9].  This  method,  however,  is  im¬ 
practical  because  of  limitations  in  present-day  CGH  fabrication  tech¬ 
nology.  Recording  devices  that  are  used  to  fabricate  CGHs  are  lim¬ 
ited  in  resolution  to  spatial  frequencies  that  are  far  below  those 
required  for  the  combiner.  CGHs ,  \noreover ,  are  not  formed  as  volume 
holograms,  and  they  generally  prockice  unwanted  diffracted  orders  in 
addition  to  the  desired  diffracted  order. 


By  combining  a  CGH  with  conventional  optical  elements,  we  can 
make  use  of  its  capability  for  forming  arbitrary  phase  variations 
that  might  otherwise  be  unrealizable.  For  this  method,  the  desired 
construction  wavefront  at  the  combiner  is  formed  by  an  optical  system 
in  which  conventional  elements  provide  the  basic  shape  and  the  CGH 
adds  small  but  irregular  phase  perturbations  to  the  wavefront.  In 
addition,  the  conventional  elements  are  arranged  so  that  there  is  a 
plane  where  a  spatial  filter,  or  aperture,  can  be  located  to  block 
unwanted  orders  diffracted  by  the  CGH. 


The  task,  then,  of  designing  a  construction  beam  is  as  follows. 
Starting  with  a  spherical  wavefront  diverging  from  a  point  source, 
an  optical  system  is  required  that  will  provide  the  desired  phase 
variation  at  the  hologram  surface,  and  at  an  intermediate  location 
form  a  real  image  of  the  point  source  where  a  spatial  filter  can  be 
located.  Figure  31  shows  the  form  that  such  an  optical  system  might 
have. 

Since  we  are  given  the  wavefront  description  at  the  hologram,  we 
can  analyze  the  construction  beam  in  the  reverse  direction,  tracing 
rays  away  from  the  hologram  toward  the  point  source.  For  a  system 
having  the  form  shown  in  Figure  31,  the  objective  is  to  select  con¬ 
ventional  optical  elements  that  focus  the  wavefront  to  the  interme¬ 
diate  point  focus  at  the  spatial  filter  plane,  and  then  to  a  second 
point  focus  that  will  be  the  location  of  the  construction  beam  point 
source.  The  quality  of  the  intermediate  point  focus  need  only  be 
sufficient  to  permit  separation  of  the  orders  diffracted  by  the  CGH. 
The  point  focus  at  the  point  source  location  also  does  not  need  to 
be  of  high  quality  as  the  final  correction  will  be  provided  by  the 
CGH.  The  point  focus  must  be  good  enough,  however,  so  that  the 
amount  of  correction  required  by  the  CGH  is  within  its  information 
capacity. 

2.  CGH  BANDWIDTH 

The  major  limitation  of  the  CGH  is  the  limited  bandwidth  of  the 
information  it  can  record.  For  a  Burch  type  CGH  [10]  recorded  at  a 
resolution  (or  spot  size)  of  ax,  the  maximum  bandwidth  af  of  the  in¬ 
formation  that  can  be  recorded  is 


For  the  Optronics  Model  1600  film  recorder  at  ERIM,  ax  =  50  pm  and 
af  =  10  cycles/mm.  For  a  recorder  that  has  a  5  pm  spot  size,  af  = 
100  cycles/mm,  which  is  a  reasonable  upper  limit  for  current 
recorders. 
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Referring  to  Figure  31,  the  CGH  converts  a  spherical  wavefront 
into  the  wavefront  required  by  the  construction  beam  design.  The 
required  CGH  bandwidth  can  be  found  from  an  analysis  of  the  diffrac¬ 
tion  required  by  the  hologram.  We  can  express  the  grating  equation 
in  the  form 


j-  Od  -  U  =  x  (48a) 

x 

-I—  (m.  -  m^)  =  x  (48b) 

y 


where  f  and  f  are  the  grating  spatial  frequencies  in  the  x  and 
x  y 

y  directions,  respectively,  1  and  m  are  ray  direction  cosines,  and 
the  subscripts  i  and  d  represent  the  incident  and  diffracted  rays, 
respectively.  Solving  for  f  and  f  ,  we  have 

*  y 


f 


X 


(49a) 


f 


y 


(49b) 


By  tracing  rays  from  the  hologram  to  the  CGH  plane,  we  can  find  the 
direction  cosines  of  the  rays  after  diffraction,  and  knowing  the  lo¬ 
cation  of  the  point  source  gives  the  direction  cosines  of  the  inci¬ 
dent  rays.  If  we  assume  that  the  CGH  lies  in  the  xy-plane,  and  that 
the  point  source  is  on  the  z-axis  a  distance  zQ  from  the  CGH,  the 
direction  cosines  of  the  incident  rays  at  a  point  (x,y)  on  the  CGH 
are 


where  R  = 
written  as 


(X2  +  y2  +zg)^2.  Equations 


(49)  can  then  be 


*  _ 
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(50a) 


fx(«.y)  -{[yx.y)  -  J] 

fy(x,y)  =  j  [md(x,y)  -  £  ]  (50b) 

The  required  bandwidths  af  and  af  of  the  CGH  can  be  determined 

x  y 

by  searching  over  the  CGH  to  find  the  maximum  and  minimum  spatial 
frequencies  required,  with 


af 

y 


^y^max  ~  ^y^min 


In  a  previous  investigation  reported  in  Reference  1,  we  used  a 
CGH  that  had  maximum  bandwidths  in  both  x  and  y  directions  of  less 
than  1  cycle/mm.  Analysis  of  a  construction  beam  similar  to  the 
final  combiner  reference  beam  indicated  that  without  the  auxiliary 
conventional  optics,  the  required  CGH  bandwidths  would  be  afx  =  70 
cycles/mm  and  afy  =  270  cycles/mm.  The  addition  of  conventional 
optics  to  provide  basic  wavefront  shaping  reduced  the  required  CGH 
bandwidths  to  30  and  20  cycles/mm,  respecti vely;  these  fall  within 
the  capabilities  of  existing  recording  devices. 


93 


* 


SECTION  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  design  investigation  of  a  low-distortion  HUD  optical  system 
demonstrated  the  usefulness  of  aspheric  holograms  that  have  analyti¬ 
cally  defined  construction  wavefronts.  These  holograms  not  only 
offer  a  great  deal  of  design  flexibility,  but  also  permit  the  de¬ 
signer  to  concentrate  on  the  hologram  itself,  postponing  the  design 
of  the  construction  beams  until  the  hologram  design  is  complete. 
Previous  design  methods  have  required  complete  specification  of  both 
construction  beams,  including  all  required  optical  components,  prior 
to  the  hologram  ray  trace.  A  key  tj  the  application  of  analytically 
defined  construction  wavefronts  is  the  use  of  computer  generated 
holograms,  which  greatly  increases  the  flexibility  of  the  construc¬ 
tion  wavefront  design,  permitting  the  specification  of  wavefronts 
that  would  be  difficult  or  impossible  to  realize  with  conventional 
optics  alone. 

From  our  design  investigation  of  a  low-distortion  HUD  optical 
system,  we  conclude  that  distortion  correction  to  the  goal  of  0.5 
percent  is  feasible  using  the  approach  of  the  aspheric  holographic 
combiner.  Although  our  design  did  not  meet  the  goal  over  the  entire 
field-of-view,  we  believe  that  with  additional  effort,  we  would  have 
continued  to  obtain  improved  performance.  More  work  was  also  re¬ 
quired  on  the  design  to  reduce  the  errors,  especially  the  parallax 
errors.  Some  modifications  to  the  merit  function  structure  would  be 
useful  to  provide  for  direct  parallax  computations,  so  that  the  par¬ 
allax  errors  would  receive  greater  attention  during  the  optimization 
process.  In  order  to  concentrate  on  evaluating  the  capability  of 
our  approach  to  correct  distortion,  we  limited  the  design  optimiza¬ 
tion  to  the  instantaneous  field-of-view  (12’v  x  20°H).  Further  op¬ 
timization  of  the  design  should  be  expanded  to  cover  the  total  field- 
of-view  (20<>V  x  20°H).  We  found,  on  the  other  hand,  that  once  the 
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initial  aspheric  combiner  design  had  been  established,  further  opti¬ 
mization  of  the  combiner  design  alone  was  not  fruitful.  The  most 
meaningful  results  were  obtained  from  optimization  of  the  system  as 
a  whole,  based  on  complete  system  performance. 

For  the  design  we  have  reported  on  in  this  report,  further  work 
is  recommended  along  the  lines  described  above.  We  also  recommend 
continued  investigation  of  aspheric  holograms  at  a  more  basic  level. 
The  investigation  of  holographic  imagery  from  the  point-of-view  of 
minimizing  the  wavefront  error  as  described  in  Section  II  warrants 
continued  study.  The  investigation  of  other  methods,  such  as  Zernike 
polynomials,  for  describing  the  construction  wavefronts  may  result 
in  significant  improvements  in  hologram  flexibility  and  performance. 
Of  interest  specific  to  the  HUD  application  would  be  an  investigation 
of  coma  correction  in  the  combiner  by  means  of  improved  construction 
wavefront  design,  as  coma  is  the  dominant  aberration  after  correction 
of  astigmatism,  and  furthermore,  coma  correction  in  the  relay  lens 
accounts  for  much  of  the  decentering  and  tilting  of  the  relay  lens 
elements.  It  would  also  be  of  interest  to  investigate  the  introduc¬ 
tion  of  one  or  more  additional  aspheric  holographic  elements  into 
the  optical  system  with  the  objective  of  simplifying  the  conventional 
optics  portion  of  the  relay  lens. 
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